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Abstract 

We perform a systematical and analytical study of lepton mixing which can be 
derived from the subgroups of SU (3) under the assumption that neutrinos are Dirac 
particles. We find that type D groups can predict lepton mixing patterns compatible 
with the experimental data at 3a level. The lepton mixing matrix turns out to be of 
the trimaximal form, and the Dirac CP violating phase is trivial. Moreover, we extend 
the flavor symmetry to the quark sector. The Cabibbo mixing between the first two 
generations of quarks can be generated by type D groups. Since all the finite subgroups 
of U( 3) that are not the subgroups of 517(3) have not been classified, an exhaustive 
scan over all finite discrete groups up to order 2000 is performed with the help of the 
computer algebra system GAP. We find that only 90 (10) groups for Dirac (Majorana) 
neutrinos can generate the lepton mixing angles in the experimentally preferred ranges. 
The lepton mixing matrix is still the trimaximal pattern and the Dirac CP phase 
remains trivial. The smallest groups that lead to viable mixing angles are [162,10], 
[162,12] and [162,14], For quark flavor mixing, the correct order of magnitude of 
the Cabibbo-Kobayashi-Maskawa matrix elements can not be generated. Only the 
Cabibbo mixing is allowed even if we impose very loose constraints 0.1 < | V^ s | <0.3 
and \ V u b\ < |D c b| < | V us \. The group A(6x7 2 ) can predict a Cabibbo angle 8 q = 7r/14 in 
good agreement with the best fit value. The observed Cabibbo mixing angle can easily 
be accommodated if the first two left-handed quark fields are assigned to a doublet. 
The groups that can give rise to both phenomenologically viable lepton mixing angles 
and acceptable Cabibbo angles are discussed, and the groups A(6 x 9 2 ), [648,259], 
[648, 260], [648, 266], and A(6 x 14 2 ) are especially promising in the case of the triplet 
assignment for both quark and lepton sectors. The three groups [496,19], [496,21], 
and [496, 23] are interesting candidates for flavor symmetry if the left-handed lepton 
and quark fields are assigned to triplet and doublet plus singlet, respectively. 
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1 Introduction 


It has been firmly established that flavor mixing exists in both quark and lepton sec¬ 
tors, and the quark and lepton mixing structures are described by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix and the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, re¬ 
spectively [ ]. It is well known that the CKM matrix is nearly diagonal and its off-diagonal 
elements are rather small. On the other hand, each entry of the PMNS matrix is of order 
one except the (13) element, which is related with the reactor mixing angle 9 13 [2-4], Fur¬ 
thermore, the CP violating phase in the quark sector has been precisely measured, while the 
value of the leptonic CP phase is not known so far. Obviously the mixing patterns in quark 
and lepton sectors are quite different. This is the so-called flavor puzzle of the standard 
model, yet the underlying origin is not revealed. 

There have been many attempts in the literature to explain the observed flavor mixing 
structures. The most popular approach is to invoke some flavor symmetry acting on different 
families to describe the observed patterns. In particular, the finite discrete symmetry group 
was extensively explored to understand the leptonic mixing angles (please see Refs. [5-9] for 
some recent reviews). In this setup, a discrete flavor symmetry group is broken down to 
different subgroups in the neutrino and charged lepton sectors and the mismatch between 
the two subgroups determines the PMNS matrix. Conversely, the flavor group can be re¬ 
constructed from the residual symmetries of the mass matrices. This approach was first 
motivated by the distinctive tri-bimaximal mixing pattern [10]. After the reactor angle 
#13 is precisely measured to be nonvanishing by T2K [11], MINOS [12], Double-Chooz [13], 
RENO [14], and DAYA-BAY [ 3], much effort has been devoted to searching for finite flavor 
symmetry groups that can give rise to neutrino mixing angles in agreement with experi¬ 
mental data. Under the assumption of Majorana neutrinos, a comprehensive scan of SU{3) 
subgroups with order less than 512 [16] as a flavor symmetry has been performed. A more 
complete search for f/(3) subgroups of order smaller than 1536 has been undertaken [17]. 
It is amazing that only three groups A(6 x 10 2 ), (Z 18 x Z 6 ) x S 3 , and A(6 x 16 2 ) can lead 
to acceptable mixing patterns [ ]. Discrete groups of order up to 200 [18] and the “excep¬ 
tional” finite groups [19] have also been investigated with the assumption of Dirac neutrinos. 
The discrete flavor symmetry paradigm is further extended to the quark sector to predict 
quark mixing angles [18,20]. Moreover, the variant scenario of only one row or one column 
of the PMNS matrix constrained by flavor symmetry is discussed [ ]. Now there appears 
to be some consensus that at least the A(6n 2 ) groups are viable candidates for flavor sym¬ 
metry [ ]. Instead of exploiting finite non-Abelian flavor symmetries, a systematic group 
space scan of discrete Abelian flavor symmetries has been performed in Ref. [23]. 

Recent notable progress is that the complete lists of lepton mixing matrices that can be 
obtained from finite flavor symmetries have been derived for Majorana neutrinos [2 ]. It is 
found that the mixing patterns are strongly restricted, the second column (in absolute values) 
of the PMNS matrix has to be (1,1, l) T / a/ 3 in order to accommodate the experimental data, 
and the Dirac CP phase is trivial. I 11 this work, we want to explore whether new mixing 
patterns can be achieved if neutrinos are Dirac particles instead of Majorana particles. We 
analytically compute all the possible mixing patterns in the case that the flavor symmetry 
group is a finite subgroup of £77(3) or it belongs to E(3A r3 ), which is a known subgroup 
series of U( 3). In particular, the mixing patterns that can be derived from the discrete 
groups ( Z 3n x Z n ) x S 3 will be systematically studied for the first time. As the vast majority 
of finite subgroups of U( 3) are not subgroups of SU(3), and a complete classification of all 
discrete subgroups of U (3) is not available, we perform a numerical scan of mixing patterns 
derived from finite discrete groups up to order 2000 with the help of powerful computer 
algebra system GAP [25-27]. Another motivation of this work is to scrutinize whether viable 
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quark mixing angles together with leptonic mixing angles can be reproduced from the same 
discrete flavor symmetry group, and which group is the smallest one capable of achieving 
this. 

The remaining parts of this paper are organized as follows. In Sec. 2 we review the 
idea of discrete flavor symmetry from both bottom-up and top-down approaches, and the 
method of extracting the flavor mixing matrix from residual flavor symmetry is presented. 
In Sec. 3, we study the mixing patterns that can be derived from the subgroups of SU (3) 
under the assumption that neutrinos are Dirac particles. We find that the subgroups of 
type D can give rise to lepton mixing patterns that are compatible with experimental data, 
the PMNS matrix is of the trimaximal form, and the Dirac CP violating phase is trivial. 
Moreover, the flavor symmetry is extended to the quark sector, and the Cabibbo mixing can 
be generated by the type D groups. In the same fashion, the mixing patterns arising from 
the U( 3) subgroup series E(3iV 3 ) are discussed in Sec 4. However, no mixing pattern in the 
experimentally preferred 3cr range is found unless the remnant symmetries are partially or 
completely accidental. Since there are a lot of finite groups that are subgroups of 17(3) but 
not subgroups of 517(3), in Sec. 5 we perform an exhaustive scan of finite discrete groups 
with order less than 2000 with the help of the powerful group software GAP. The possible 
residual symmetries and the corresponding mixing matrices as well as mixing parameters 
for each group are available at the Website [ ]. We find that the lepton flavor mixing in 

accordance with experimental data is still of the trimaximal pattern, and only the Cabibbo 
mixing is obtained even if we impose very loose constraints on the CKM matrix elements 
0.1 < \V US \ < 0.3 and \V u b\ < \V c b\ < \V US \. We summarize our results and conclude this 
paper in Sec. 6. Finally the procedure of how to convert a nonunitary representation into a 
unitary representation is presented in Appendix A. 

2 Basic framework 

In this section, we briefly review the basic idea of flavor symmetry and how the mixing 
matrices can be predicted from the structure of flavor symmetry group Qj? and the assumed 
symmetry breaking pattern. In the bottom-up approach, the mass terms of the Lagrangian 
for quarks and charged leptons is of the form 

Cfnass — —UrttlijUl — D R mr>D L — IrwliIl + H.c ., (2.1) 

where U L = (u L ,c L ,t L ) T , D L = ( d L ,s Lj b L ) T , l L = (e L , /i L , r L ) T J and Ur = ( u R ,c R ,tR) T , 
Dr = (<7r, Sr, bR) T , Ir = (eR, hr,tr) t represent the three generations of left-handed and 
right-handed up type quark, down type quark, and charged lepton fields. If neutrinos are 
Dirac particles, the mass terms are 

£ V mass = -VRinvVL + H.C. , (2.2) 

where u L = (u eL , v rL ) T and u R = (u eR , u^r, v tR ) t are the left- and right-handed neutrino 
fields. For Majorana neutrinos, it is 

C s = \vlC- l m u v L + H.c., (2.3) 

where C is the charge-conjugation matrix. Both quark and lepton mass matrices are diago¬ 
nalized by 


V^m^muVu = diag (m 2 u , m 2 c ,m 2 t ) , 
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(2.4) 


V ] D rrJ D m D V D = diag (m 2 d , m 2 , m 2 b ) , 
V^mjmiVi = diag {m 2 , m 2 ^, ml) , 


and 


V l }mlm u V l/ = diag , for Dirac neutrinos, 

V^m v V v = diag m 3 ), for Majorana neutrinos. (2.5) 


The CKM and PMNS matrices are defined as 

Vckm = V^Vd, Dpmns = V^I4. (2.6) 

One can straightforwardly check that the Hermitian combinations m'jjmu, m^mn, and m\mi 
are invariant under the unitary transformations, 


Ul^-TjjUl, D l T d D l , l L ^-Til L , (2.7) 

with 

T c = V c diag (e ia % e i/3c , e^) V? , (2.8) 

where the subscript c G {U,D,l}, a c , (3 C , and y c are arbitrary real parameters, i.e., 

Tlm\m c T c = m\m c . (2.9) 

In the language of group theory, each of the quark and charged lepton mass terms has a full 
17(1) x 7(1) x U( 1) symmetry. In the same fashion, for the Dirac neutrino, we have 

vl ->■ T v v l , mlm u -x Tlmlm u T v = mj^m„ , ( 2 -10) 

with 

T v = Kdiag (e ia \ Vj , (2.11) 

where f3 u , and are real. In the case of the Majorana neutrino, the symmetries of the 
neutrino mass term in Eq. (2.3) are strongly constrained. It has been established that the 
symmetry transformation of m v is [10] 


->■ T v u l , m u -x T^m u T v = m v , (2.12) 

with 

T v = 14diag (±1, ±1, ±1) Vj . (2.13) 

Note that if the lightest neutrino is massless, the corresponding diagonal entry “±1” could 
be an arbitrary phase factor. The eight possible values of T u in Eq. (2.13) constitute the 
symmetry group Z 2 x Z 2 x Z 2 . However, only Z 2 x Z 2 is effective since the third Z 2 just 
contributes an overall “—1” factor that can be absorbed by redefining the neutrino fields. 
In other words, the effective symmetry of the neutrino mass term is the Z 2 x Z 2 Klein four 
group for the Majorana neutrino. Inversely, given the symmetry transformations that leave 
the quark and lepton mass terms invariant, one is able to reconstruct the explicit form of 
their mass matrices. 

In the top-down approach, the above mentioned symmetries of the quark and lepton 
mass terms originate from some flavor symmetry group Gj- at high energy scale. The full 
Lagrangian is invariant under Gt, and then subsequently Gt is broken down to different 
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subgroups Gu and Gi in the neutrino and charged lepton sectors, and to Gu and Gd in the 
up quark and down quark sectors, respectively. From Eqs. (2.8, 2.11, 2.13), we know 

Gc C 1/(1) X 1/(1) X 1/(1), 

G u ^ 1/(1) X 1/(1) x 1/(1), for Dirac neutrinos, 

Gu Q Z 2 x Z 2 x Z 2 , for Majorana neutrinos. (2-14) 

That is to say, the residual symmetries Gu Gu, and Gd can be any Abelian subgroups of Gt, 
and similarly Q v can be any Abelian subgroup of Gt if neutrinos are Dirac particles while Gu 
can only be effectively (or a subgroup of) a Klein four subgroup for Majorana neutrinos. In 
this work, we consider the scenario that Gt is a non-Abelian finite discrete group. As usual, 
the three generations of left-handed lepton doublets are assigned to be an irreducible faithful 
three-dimensional representation p 3 of Gt [16-19,21,22]. Analogously the three left-handed 
quark doublets are also embedded into the same three-dimensional representations [18,20]. 
For the residual symmetries Gu,i,u,d to hold, the quark and lepton mass matrices have to 
satisfy the following constraints: 

pl(9c)mlm c p 3 (gc) = m ] c m c , g c ^G c , (2.15) 

and 

pl(g u )mlm u p 3 (g u ) = m\,m v , for Dirac neutrinos, g u G Gu, 
p 3 {gu)'^uPs{gu) = mu, for Majorana neutrinos, g v G Gu ■ (2.16) 

Notice that it is sufficient to impose the invariant conditions of Eqs. (2.15, 2.16) on the 
generators of G c and G v - Furthermore, we see that Eqs. (2.15, 2.16) imply 

[p3(gc),mlm c \ = 0, [p 3 (gu),mlmu\ = 0. (2.17) 

As p 3 (g c ) and m\m c commute with each other, the unitary transformation V c that diago¬ 
nalizes the mass matrix rrv c m c also diagonalizes P 3 (g c ) up to permutations and phases of 
columns, as shown in Eqs. (2.4) and (2.8). Similarly p 3 (g v ) and m\,mu (or m„) are diago¬ 
nalized by the same matrix V v , as displayed in Eqs. (2.5, 2.11, 2.13). I 11 other words, given 
a group representation, the unitary transformations V c and V v (and eventually the CKM 
matrix Vckm and the PMNS matrix Dpmns) can be obtained by diagonalizing the represen¬ 
tation matrices of the generators of Gc and Gu, respectively. In practice, we only need to 
find the eigenvectors of p 3 (g c ) and ps{g v ) that form the column vectors of V c and V v . In 
this approach, it is not necessary to construct the explicit form of the mass matrices rru c m c 
and m\,m v (or m u ) although this can be accomplished straightforwardly. Since the order 
of both quark and lepton masses cannot be pinned down in this framework, particularly 
the neutrino mass spectrum can be either normal ordering (NO) or inverted ordering (10), 
the unitary matrices V c and V v are uniquely fixed up to permutations and phases of their 
column vectors. As a result, both C/pmns and Vckm are determined up to independent row 
and column permutations and arbitrary phase matrices multiplied from the left-handed and 
right-handed sides. Notice that if we switch the roles of the subgroups Gu and Gi, the PMNS 
matrix 7/pmns would become its Hermitian conjugate. Analogously, the exchange of Gu and 
Gd would lead to Vckm being Hermitian conjugated. 

It is remarkable that the PMNS matrix I/pmns stems from the mismatch between Gu and 
Gi and their relative embedding into Gt, and Vckm from the misalignment of Gu and Gd in 
this framework. Physical results only depend on the structure of flavor symmetry group Gt 
and the assumed symmetry breaking patterns, and they are independent of how the required 
residual symmetries are dynamically realized. 
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3 Flavor mixing achievable from subgroups of SU( 3) 

The finite subgroups of 517(3) have been systematically classified nearly 100 years ago [29]. 
They are further studied in Refs. [7, 30-31 ] in recent years, and the complete list can be 
found, for example, in [7,33]. Now it has been firmly established that all discrete subgroups 
of 57/(3) can be divided into five categories: type A, type B, type C, type D, and type E. 
In the following, we shall investigate the viability of each type of these groups as a flavor 
symmetry for Dirac neutrinos and further study the predictions for the mixing matrix. The 
groups of type A are Abelian groups, and they are isomorphic to Z m x Z n , where n is a 
divisor of m. Obviously type A subgroups are not suitable as flavor symmetry groups due 
to the Abelian characterization. Groups of type B are finite subgroups of 17(2), and up to 
basis transformations, such 577(3) matrices have the form 

( d C °m) w,th MeC7(2) ' ^ 

This means that the three-dimensional representations are reducible. As a consequence, 
at least four entries of the mixing matrix would be vanishing such that the reactor mixing 
angle 6\ 3 would be zero. Therefore type B subgroups can-not give phenomenologically viable 
lepton mixing patterns. Type E subgroups are the so-called “exceptional” finite groups of 
517(3), and they include £(60), £(60) x Z 3 , £(168), £(168) x Z 3 , £(36 x 3), £(72 x 3), 
£(216 x 3), and £(360 x 3) [7,3 ]. The mixing patterns that can be derived from these 
exceptional groups have been comprehensively studied in Ref. [19]. The results of [19] are 
confirmed by our systematical numerical analysis for the groups up to order 2000 in the next 
section. Interestingly, the groups of type C comprise infinitely many series, two of which are 
the well-known A(3n 2 ) and T m series. Nevertheless, all groups of type D are described by just 
two infinite series where one is the well-known series A(6n 2 ). Notice that A(6 n 2 ) as flavor 
symmetry groups have been widely discussed [22,34-37], and they remain viable candidates 
for flavor symmetry [17,18,38]. In this section, we shall present a complete analytical study 
of mixing patterns obtained from type C and type D groups. 

As shown in Refs. [33,39], the structures of groups of type C and D are given by 

type C : Q= (Z m x Z n ) x Z 3 , 
type D : Q= (Z m x Z n ) x S 3 , 

where S 3 is isomorphic to Z 3 x Z 2 , and n is a divisor of m. It is worth mentioning that the 
above group structure crucially depends on the definition of the semidirect product. A new 
group could be obtained if the definition of the semidirect product is changed. We really find 
that some U( 3) subgroups apparently have the same structure as the type D groups while 
they are, in fact, different because of different semidirect structures. The groups of type C 
can be conveniently expressed in terms of three generators a, c, and d, while type D groups 
involve one more generator b, and the explicit forms of these generators are 

/0 10\ /-I 0 0 \ (rj 0 0 \ /I 0 0\ 

a =j00l), 6 =( 0 0 —1 ) , c = I 0 77 * 0 , d = 0 ?R r 0 , (3.3) 

\10 0/ \ 0 —1 0 / \0 0 \0 0 rf) 

where rj = e 2m / m . In fact, the matrices in Eq. (3.3) form a faithful irreducible representation 
of the type C and type D groups. One can check that the following multiplication rules are 
satisfied: 

a 3 =b 2 = ( ab ) 2 = c m = d n = 1, cd = dc, 
aca -1 = c k d l , ada~ l = c~ r d~^ k+1 \ (3.4) 

bcb~ l = cd 1 , bdb~ l = d ~ l , 
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Group series 

Constraints on parameters 

ri ( fc ) 
rn,n 

3r'n,n ^3 * r'n,n 

'~'9r'n' ,3n' 

m = rn, 3 \ r 

77i = 3r'n, r = 3 r', 3 \r', 3 \n 
m = 9r'n', n = 3n', r = 3r', 3 \ r' 

D ( °l = A (6 n 2 ) 
D^ n = Z 3 x A(6n 2 ) 

^9n',3n' 

m = n, r = 1 
in = 3 n, r = 3, 3 \ n 
m = 9n', n = 3n', r = 3 


Table 1: Summary of the groups of type C and type D. The notation \ denotes “does not divide.” For the 
type C groups, the parameters r and r' can be any product of primes of the form 6^ + 1 with non-negative 
integer z. The value of k is fixed by Eq. (3.5). 


where r = m/n, and the last two lines define the semidirect product. For groups of both 
type C and type D, the integers k, /, and r fulfill the relation: 

1 + k + k 2 = Ir . (3.5) 

For groups of type D, the parameters k, l', and r satisfy 

1 + 2 k = Hr . (3.6) 

As we can see, groups of type C and type D are completely fixed by the parameters m, n and 
k according to the multiplication rules in Eq. (3.4). Therefore, we can denote these groups 
by the symbols Cm} n and Dm, respectively. 

The group structures of type C and type D groups have been studied in Ref. [33]. The 
complete classification of these groups is summarized in Table 1. Notice that there is an 
infinite series Crn,n for every pair (r, k). Therefore infinitely many type C series are obtained 
since there are infinitely many admissible r. In Table 2 we have listed all allowed values of 
r with r < 100 and the corresponding solution(s) of k. Some well-known type C groups are 
included in table 2. For example, if we fix r = 1 and k = 0, we shall get the dihedral-like 
groups 

Cg = A(3n 2 ) ^ (Z n x Z n ) x Z 3 . (3.7) 

If we set n — 1, r = m is the product of primes of the form 6^ + 1 with z being a positive 
integer, and then the so-called T m group is obtained, 

Cffi = T m = Z m x Z 3 , (3.8) 

where m and k are related by 

1 + k + k 2 = 0 (mod m ). (3.9) 

From Table 2, we see that these groups exist only for specific integers m = 7, 13, 19, 31, 37, 
43, 49, 61, 67, 73, 79, 91, 97, ..., and there are two solutions of k = 9 and 16 for m = 91. 
Consequently there are two nonisomorphic groups T 91 and T 91 . 

For type D groups, the values of the parameters r and k are strongly constrained. It is 
found that only two cases are admissible: r = 1, k = 0, m = n or r = 3, k = 1, m = 3n. 
Therefore all type D groups can be classified into two series Dn.i and D^ n . Notice that 
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r 

k 

r 

k 

1 

0 

7 x 7 = 49 

18 

3 

1 

3 x 19 = 57 

7 

7 

2 

61 

13 

13 

3 

67 

29 

19 

7 

73 

8 

3 x 7 = 21 

4 

79 

23 

31 

5 

7 x 13 = 91 

9,16 

37 

10 

3 x 31 = 93 

25 

3 x 13 = 39 

16 

97 

35 

43 
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Table 2: List of all allowed r < 100 and the corresponding solution(s) k of Eq. (3.5). 


Dn} n is exactly the well-known A(6n 2 ) group [10], and the chosen generators a', b' , c 1 , and 
d! of Ref. [40] are related to the present ones by 

a = a, b' = ab, c = cd , d' = d -1 . (3.10) 

Furthermore, D^ n is isomorphic to Z 3 x A(6n 2 ) if n is not divisible by 3 [ ]. As a result, 

all type D groups can be genuinely described by A(6n 2 ) and D^j, 3n , = (Z 9n r x Z 3 n i) x S 3 . 

3.1 Remnant symmetries and their diagonalization 

In this section, we shall discuss the possible remnant symmetry groups and the corre¬ 
sponding unitary diagonalization matrices, if the flavor symmetry group Qjr is a group of type 
C or type D. As shown in Sec. 2, both residual symmetries Q u and Qi are Abelian subgroups 
of Qjr if neutrinos are Dirac particles. The same is true for Qu,d in the quark sector. First we 
consider the situation that the remnant symmetry group Q r with r = u, l, U, D is generated 
by a single element. If the generator of Q r is represented by a matrix with nondegenerate 
eigenvalues, the diagonalization matrix V r of Q r is uniquely determined by its eigenvectors 
up to permutations and phases. On the other hand, if two of the three eigenvalues are 
degenerate, only one generation can be distinguished from the other two generations and 
consequently only one column of V r can be fixed by the postulated residual symmetry. Such 
a column will be called “mixing vector” hereafter. The generator of Q r can be written as 
c s d*, 6c s d*, ac s d*, a 2 c s d t , a6c s d*, or a 2 6c s d* with s — 0,1,..., m — 1, t — 0,1,..., n -=• 1. The 
corresponding unitary transformation V r can be straightforwardly calculated as follows: 

• Q r = ( c s d l ) 

The representation matrix of the generator c s d* is given by 

(r) s 0 0 \ 

P 3 {c s d t ) = 0 T] ks ~ rt 0 

\ 0 0 ^-(fc+1)* J 

Its eigenvalues are 

s ks-rt rt-(k+l)s 

'I 1 'I J 'I 


(3.11) 

(3.12) 
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The unitary matrix V r , which diagonalizes P 3 (c s d*), is of the form 


M : V r 



(3.13) 


Nondegeneracy of the eigenvalues imposes the following constraints on the parameters s 
and t\ 


(2k + l)s — 2 rt ^ 0 (mod m), 

(k + 2)s — rt ^ 0 (mod m), 

(k — l)s — rt ^ 0 (mod m ). (3-14) 

In the case of two degenerate eigenvalues, only one column of V r can be fixed. This 
column vector would be denoted by v r henceforth. To be more precise, we have 

V i : v r = (1, 0, Of for (2k + l)s — 2 rt = 0 (mod in), (k + 2)s — rt ^ 0 (mod m), 

T >2 : v r = (0,1, Of for (k + 2)s — rt = 0 (mod m), (k — l)s — rt ^ 0 (mod m), 

V 3 : v r = (0, 0, if for (k — l)s — rt = 0 (mod m), (k + 2)s — rt ^ 0 (mod m) . (3.15) 


Gr = (■ be s d t ) 


-r] s 0 0 

p 3 (bc s d *) = | 0 0 - v rt~(k+i)s | ) 

0 -rjks-rt 0 


with the eigenvalues 
The diagonalization matrix V r is given by 


- V 3/2 , rj 3/2 , -rf . 


0 0 y/2' 

J\f 2 : V r = f= | rj rt -( k +D s -rft-( k +h)s q 


1 

ft 


1 1 0 
To avoid degeneracy among the eigenvalues, we obtain the constraint 

3s 7 ^ 0 (mod m). 

If two of the eigenvalues are degenerate, the fixed column vector is 

V 4 : f= ^0, r] rt ~^ k+ ^ s , lj for 3s = m (mod 2m), 

1 


D 5 : — ^0, —rf ( fc + 2 ) s ; lj for 3s = 0 (mod 2m). 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


• G r = ( ac s d b ) 

The order of the generator cur'd* is always three; therefore residual symmetry Gr = 


is a Z 3 subgroup with 


0 r] ks ~ rt 0 

p 3 (ac s cf = [ 0 0 rf-tfc+rp 

r/ s 0 0 


(ac s d*) 


(3.21) 



which is diagonalized by 


1 / TJ s UJ7] s L0 2 7] s \ 

A /3 : V r = — ^(fc+lja w 2 ?? rt-(fc+l), ur] rt-(k+l)s # (3.22) 

V 3 Y 1 1 1 / 

where cu = e 2 * 71 "/ 3 . The eigenvalues are 1, cu, and or for any values of s and t. 

• G r = (a 2 c s $) 

In this case, Q r = (a 2 c s c/*) is also a Z 3 subgroup for any values of s and t, 

/ 0 0 77 rt -( fe+1 ) s \ 

Pz{a 2 c s d t ) = 77 s 0 0 , (3.23) 

\ 0 T) ks ~ rt 0 / 

whose eigenvalues are 1, u, and to 2 . Since the eigenvalues are not degenerate, the unitary 
matrix V r is completely fixed 


V ^7S 


' ifjTt— (fc+l)a ^^rt-(fc+l)s ^2^rt-(fc+l)s 


~rt—ks 


^j2 ^jrt—ks 


UT] 


rt—ks 


(3.24) 


Note that this V r can be obtained from the one in Eq. (3.22) via the replacement s —> 
(k + l)s — rt, t —* s(l + k + /c 2 )/r — kt. This fact is because of the relation (a 2 c s e / f ) 2 = 

ac (fc+l)s—rt^s(l+fc+fc 2 )/r—fct 

• I?,- = (ab^d*) 


0 0 

ps(abc s d t ) = ( 0 -r/ fcs - rf 

- 77 s 0 


_ r jrt—(k+ l)s ' 

0 
0 


which is diagonalized by the unitary transformation 

0 ^fhi-(fe+2)s) _ , n \(rt-{k+2)s) 


7/s 


A4 : K - — | y/2 


0 


0 

1 


0 

1 


with the eigenvalues 


V 


rt—ks ) 


_ nks-rt _ § (rt-ks) 

7 ) 7 > 

The nondegeneracy of the solutions requires 

3 (rt — ks ) 7 ^ 0 (mod m). 

The degenerate solutions are given by 

V G : v r = -^= 0 , f or 3(rf — /us) = m (mod 2m), 

V 7 :v r = ^= ( 


V2 V 


, 0,1 j for 3(rt — /us) = 0 (mod 2m). 


(3.25) 


(3.26) 


(3.27) 

(3.28) 


(3.29) 
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Q r = { K a 2 bc s d t ) 


0 _■ rjks-rt 0 

p s (a 2 bc s d lt ) = [ —if 0 0 

0 o _^-(A;+i)s 

The unitary diagonalization matrix V r is determined to be 


1 / 0 d 

A 5 : V r = — 7 = I 0 1 

V2 ' V 2 0 


5 ((k—l)s—rt) 


V 


\{(k-l)s-rt) 

l 


,i((fc+l)s-rf) i((fc+l)s-rt) 


The corresponding eigenvalues are 

- r/ ri -( fc+1)s , -772' 

with the nondegenerate condition 

3 ((k + l)s — rt ) 7 ^ 0 (mod m). 

The partial degeneracy of the eigenvalues leads to 

V 8 : v r = ~^= ^r/ 2 ^ k ~ 1)s ~ rt \ 1, oj for 3((fc + 1 )s — rt) = m (mod 2m), 

£>9 : v r = ~^= 1 , 0 ^ for 3 {{k + l)s — rt) = 0 (mod 2 m). 


(3.30) 


(3.31) 


(3.32) 


(3.33) 


(3.34) 


As shown above, only one row of V r can be fixed for the cases with degenerate eigenvalues 
such as D\, ..., Vg. To resolve this degeneracy and fully determine the unitary matrix V r , we 
extend the remnant symmetry group to be the product of cyclic groups, e.g., Q r = G 1 x Q 2 . 
Here we follow the minimality principle in which either Gi or G 2 alone is insufficient to 
distinguish the three generations of fermions. 1 That is to say, the generators of both Gi and 
G 2 are represented by matrices that have two degenerate eigenvalues; in such instances, one 
of Gi and G 2 alone determines only a column in V r and the third column can be fixed by the 
unitary condition. As a consequence, neither Gi nor G 2 can be the group generated by 
or a 2 c s db Moreover, by definition, the generators of Gi and G 2 should be commutable with 
each other. All admissible combinations of Gi and G 2 satisfying the above requirements are 
summarized in Table 3. 

For the cases of the degeneracy can be eliminated by taking G 1 = {c s d *) and G2 = 

( c s (f ), where s, t and s', t' should fulfill any two different conditions shown in Eq. (3.15). 
Obviously the unitary transformation V r , which diagonalizes both Gi and G2, would be a 
trivial unit matrix up to permutations and phases of rows and columns. Then we turn 
to the scenario of T> 4 ^. For G\ = (bc s d l ) with 3s = 12m (mod 2 m) and l 2 = 0 ,1, only 

one column of V r would be determined to be ^0, — (— l) l2 rj rt ~( k+ ^ s , 1^ . We can enlarge 

the residual symmetry group to Gr — Gi X G 2 with G 2 = (c^d 1 ') or G 2 = (bc^d 1 '). I 11 
case of G 2 = {c s d l ), the commutability between G 1 and G 2 leads to the condition (2 k + 
l)s' — 2 rt' = 0 (mod m). Consequently another column fixed by G 2 is (1,0, 0) T , which 
is exactly T> 1 . Therefore the unitary transformation V r is of the form of J\f 2 while the 

1 If £1 (or Q 2 ) is sufficient to distinguish among the generations, the presence of the group Q 2 (or Qf 
would not add any information as far as lepton mixing is concerned, and the unitary matrix V r would be of 
the form of Eqs. (3.13), (3.18), (3.22), (3.26), or (3.31). 
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Gi 


G 2 


Constraints on group parameters 


Form of V r 


(c s d') 


(c s 'd 1 ') 


or 


or 


(2k + l)s — 2rt = 0 (mod m) 
( k + 2) s' — rt' = 0 (mod m) 
(2k + l)s — 2 rt = 0 (mod m) 
(k — l)s' — rt' = 0 (mod m) 
(k + 2)s — rt = 0 (mod m) 
(k — l)s' — rt' — 0 (mod m) 
(s s' , t t') 


A/i 


(he'd!) 


’d?) 


(2k + l)s' — 2 rt' = 0 (mod m) 
3s = 0 (mod m) 


(bc s 'd?) 


(2k + l)(s — s') — 2r(t — t') = lim (mod 2m) 
3s = l 2 m (mod 2m), 3 s' = l^m (mod 2m) 


A fo 


(i abc s d t ) 


’d?) 


(k + 2)s' — rt' = 0 (mod m) 
3 (rt — ks) = 0 (mod m) 


(abc s 'd t ') 


(k + 2)(s — s') —r(t — t') = 1 1 m (mod 2m) 
3 (rt — ks) = l^m (mod 2m) 

3 (rt' — ks') = l^m (mod 2m) 


Ma 


■Gx^d 1 ) 


(c s 'd 1 ') 
(a 2 bc s '($') 


(k — 1)s' — rt' = 0 (mod m) 

3 ((k + l)s — rt) = 0 (mod m) 

(k — l)(s — s') — r(t — t') = l\m (mod 2m) 
3 ((k + l)s — rt) = l 2 m (mod 2m) 

3 ((k + l)s' — rt') = l 3 m (mod 2m) 


A/5 


Table 3: The possible extension of the remnant flavor symmetry Q r = Q\ x Q 2 and the corresponding 
unitary transformation V r . Here we confine ourselves to the minimal residual symmetry Q r that allows the 
distinction of three generations. In other words, Q\ or Q 2 alone fixes only a column in V r , and the column 
vectors associated with Qi or Q 2 should be different. As a result, we have the parameters ^ 1 , 2,3 = 0,1 and 
li + h +13 = 1 , 3 . 

constraint 3s 7 ^ 0 (mod m) is removed. For the second possibility G 2 = (bc s 'd 1 ') with 
3s' = l 3 m (mod 2m) and l 3 = 0,1, the product Q 1 x Q 2 is well defined if and only if we have 
(2k + l)(s — s') — 2 r(t — t') = l\m (mod 2m) with l\ = 0,1. Thus the mixing vector dictated 

by this G ‘2 becomes ^ ^0, — (— \) l 3'q rt '-d t +\) s ' ^ \ j = _L ^0, — (—l)b+b^rt-(fe+|)s^ fj _ Both 

mixing vectors V 3 and T) 5 would be obtained for li + l 2 + l 3 = 1,3 such that V r is still of the 
form of M 2 - In the same fashion, the degeneracy in V e 7 and 9 can also be resolved by 
extending Q r from ( a a bc s d t ) to its product with ( c s d t ) or ( a a bc s d t ) with a = 0,1, and the 
unitary matrix V r is still given by A /4 and A/ 5 , respectively, as shown in Table 3. In short, for 
type D groups, the diagonalization matrix V r enforced by residual symmetry can take only 
five distinct forms A/ 1 ^, 3 , 4,5 in Eqs. (3.13, 3.18, 3.22, 3.26, 3.31) where both s and t are free 
of any constraints. For type C groups, the generator b is absent; therefore V r can only be of 
the form of A/"i, 3 - 
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3.2 Predictions of lepton mixing matrix for Dirac neutrino 

If neutrinos are Dirac particles, both remnant symmetries Gi and Q v can be any Abelian 
subgroups that are sufficient to distinguish the three generations of leptons. As a result, the 
charged lepton and neutrino diagonalization matrices V\ and V v can be any ones of A/^ 2 , 3 , 4,5 
when the flavor symmetry is a group of type D. The PMNS matrix is defined as 

ffi’MNS = ViV v . (3.35) 


There are altogether 5 x 5 = 25 possible combinations of V) and V Ul and they can be further 
divided into seven classes as follows. Here we would like to emphasize again that the PMNS 
matrix is determined up to independent row and column permutations and the multiplication 
of arbitrary phase matrices from the left and the right sides. I 11 the following, s and t denote 
the exponents of the generators c and d, respectively in the residual symmetry Gi, and p and 
q are the exponents of c and d in Q v . 


• w,vy = ( Vi.M) 

( 100 \ 

UpMNS =010 

\0 0 1/ 


• (^,K) = (M,A 4 ), (A/'s,A/;) 


I^PMNS 


1 

71 



(3.36) 


(3.37) 


which is the so-called democratic mixing matrix in which all elements have the same 
absolute value [II]. 


(A), V v ) — (A/i, A/i), (A/i, A/i), (A/i, A/4), (A/4, A/i), (A/), A/5), (A/5, A/)), (A/2, A/2), (A/4, A/4), (A/5, A/5) 


/! 0 0 \ 

Upmns = 0 cos 6 — sin 6 I , 

y 0 sin 6 cos 6 ) 

where 

(A/i, A/i), (A/ijA/i), (A/i, A/5), 
(A/i, A/i), (A/4, A/i) , (A/5, A/i) 

(A/2, A/2) 

(A/4, A/4) 

(A/5, A/5) 

• ( Yi,V v ) = (A/i, A/4), (A/4, A/i), (A/i, A/5), (A/5, A/i), (A/4,A/5), (A/5,A/4) 

x / \/2 0 \ 

7/pMNS = - I — 1 1 v/2 1 • 

2 \ 1 -1 V 2 / 


0 = 


4’ 


(2fc+l)(s—p) —2r(t—g) 


7T, 


2m 

(fc+2)(s-p)-r(t-g). 


2m 

(1 -k)(s-p)+r(t-q) 
2m 




(3.38) 


(3.39) 


(3.40) 


Obviously all the above four possible forms of ZTpmns hr Eqs. (3.36, 3.37, 3.38, 3.40) are 
not compatible with the present experimental data [4]. 
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• (Vl,V u ) — (A /2 5 A/ 3 ) , (A/ 4 , A/ 3 ), (A/ 5 , A/ 3 ) 


C/pMNS 


1 

71 


/ \/2 cos # 
1 

\ — \/2sin 9 


— \/ 2 cos(# — |) 
1 

V^sin (0 - f) 


— 72 cos (9 + |) 
1 

\/ 2 sin (0 + |) 


(3.41) 


where 


C - 2 (t+l), + ( 2 W-l),- 2 K.-») T , 

«= < (A4.A4) 

l (JV 6 ,JV 3 ) 


(3.42) 


We see that one row of the PMNS matrix is (1,1, l)/73 in this case. The lepton mixing 
angles can easily be read out as 


sin 2 0 i 3 = 4 C os 2 ((9 + |), 


sin t /12 


9 12 = 


2 cos 2 (# — |) 

1 + 2 sin 2 (# + f) 


sin 


#23 = 


1 


1 + 2 sin 2 (<9 + f) 


( 3 . 43 ) 


Eliminating the parameter 6 *, correlations among the three mixing angles follow imme¬ 
diately: 


sin 2 0 i 2 = - ± - tan 6 >i 3-\/2 - tan 2 # 13 , 


3 sin 2 6*23 cos 2 # i3 = 1 


(3.44) 


Inputting the 3 cr range 0.0188 < sin 2 6*13 < 0.0251 [1], we obtain 0.387 < sin 2 9 12 < 0.403 
(or 0.597 < sin 2 6*12 < 0.613 for the plus sign) and 0.340 < sin 3 623 < 0.342, which are 
not favored by the present data [ ]. Alternatively, we can exchange the second and third 
rows of the PMNS matrix in Eq. (3.41), the atmospheric mixing angle changes from 623 
to 7t/2 — #23, and consequently we have 0.658 < sin 3 #23 A 0.660. Although the mixing 
pattern in Eq. (3.41) fails to describe the data well, small corrections to the leading order 
results are generally expected to exist in concrete models [5,8,9]. As a consequence, this 
pattern could possibly be brought into agreement with the experimental data after higher 
order corrections are included. This mixing pattern should be interesting from the model 
building perspective. 


• (Vi, K) = (A/ 3 , A/ 3 ) 



/ 1 + 2e iS 

cos 9 

1 — 2e iS cos (9- f) 

1 - 2 e iS cos (9 + f) \ 


j 

f7pMNS = 77 

1 — 2e l5 cos 

(* + I) 

1 + 2e lS cos 9 

1 — 2 e i5 cos (# — f) 

, (3.45) 

O 1 

^ 1 — 2e* 5 cos 

( 0 -f) 

1 — 2e iS cos (9+ f) 

1 + 2 e lS cos 9 J 



where 


k(s - p) — r(t - q) 

O = -7T, 


6 = (* + 2 )(»-ri- + -g) , r (3.46) 


m m 

The possible exchange of the second and third rows of the PMNS matrix leads to another 
different mixing pattern. The lepton mixing matrices arising from other permutations 
of rows and columns can be obtained from these two patterns by redefining the param¬ 
eter 6 . We notice that the solar and the atmospheric mixing angles are closely related 
as sin 2 #12 = sin 2 #23 or sin 2 #12 = 1 — sin 2 # 23 , which are not compatible with the pre¬ 
ferred values from global fits at the 3 cr level [ ]. Furthermore, we perform a numerical 
analysis where both 9 and 6 are treated as random numbers in the interval of [— 7 r, 7 r]. 
The correlation and the allowed region of sin # 13 and sin 2 # 12 are plotted in Fig. 1. Obvi¬ 
ously the observed values of sin #13 and sin 2 #12 cannot be accommodated simultaneously. 
Therefore this mixing pattern is not viable. 
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Figure 1: The correlation between sin $13 and sin 2 9 12 for the mixing pattern in Eq. (3.45). The pink region 
represents the possible values of the mixing angles when 9 and S freely vary in the region of [— n, 7 t] . The lrr 
and 3er bounds are adopted from Ref. [ 4 ]. 


• (Vt,V v ) = (A/'s, A//), (A/'s,A4), (A/'s,A/'s) 

1 / y/2 cos 6 1 — a/2 sin 6 \ 

t/pMNS = —/t; I — V^cos (6 — |) 1 y/2 sin (0 — |) I , (3.47) 

^ \ — \/2cos (0 + 1 \/2sin (# + |) / 

where 

r ~(2*+l)iH-2(frfl) a -2r(t- g ) 7r> (^ 3> ^ 2 ) 

0 =\ - (t+ 7„t gr+2 V , W.A) (3.48) 

l ina£±|£±n2z2£) Ti (M 3 ,Ms) 

We see that one column of the PMNS matrix in Eq. (3.47) is (1,1, 1) t /a/3. This mixing 
pattern is known as the trimaximal mixing. To investigate the permutations of rows and 
columns of t/pMNS, if is appropriate to use the representation of the permutations V G S 3 
as permutation matrices, 

V M{V) = (ep ( i), e-p (2 ), ep (3) ) with ei = (1, 0 , 0 ) T , e 2 = ( 0 , 1 , 0 ) T , e 3 = ( 0 , 0 , 1 ) T . 

(3.49) 

After performing row permutation V and column permutation V', the PMNS matrix 
becomes 

f/pMNS ~t M(T)UpmnsM(V'). (3.50) 

The constant vector (1,1,1 ) T /a/ 3 should be placed in the second column to accommodate 
the experimentally measured mixing angles, and the first and the third columns can be 
interchanged. On the other hand, the three rows can be permutated in any way. For 
all the admissible permutations, we find the resulting PMNS matrix can be obtained 
from the original one in Eq. (3.47) by redefinition of the parameter 6 up to an irrelevant 
diagonal matrix multiplied from the right-hand side with entries ±1. The replacement 
rules are listed in Table 4. Without loss of generality, we can extract the predictions for 
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\ V 

(123) 

(213) 

(321) 

(132) 

(312) 

(231) 

(123) 

#^# 


#->-#-- 

3 

# -+ -# 

9^9-- 

3 


(321) 

e^e + T - 


+ 

Qi 

1 

t 

#->•-#-- 

2 

#^#+^ 

6 



Table 4: The replacement rules by which the permutated PMNS matrix can be obtained from the initial 
one in Eq. (3.47), where V and V' refer to permutation of rows and columns respectively. Here a generic 
permutation (1,2,3) —> (ni, 712 , 713 ) is simply denoted by 


the mixing parameters from the PMNS matrix in Eq. (3.47). Since each matrix element 
is real, the Dirac CP phase is trivial, i.e., sin^p = 0. The lepton mixing angles are 
given by 


sin 2 #13 = - sin 2 #, sin 2 #12 = 


1 + 2 cos 2 # ’ 


sin" #23 = 


2 sin 2 (# — |) 
1 + 2 cos 2 # 


(3.51) 


As the three mixing angles depend on a single parameter #, they are strongly correlated, 


3 sin 2 #12 cos 2 #13 = 1, 


sin 2 #23 = - ± - tan # i3 2 — tan 2 # i3 , (3.52) 


which are displayed in Fig. 2. The solar mixing angle has a lower bound given by 
sin 2 #12 > 1/3 due to the first relation in Eq. (3.52). The measurement of the reactor 
mixing angle sin 2 # 13 ~ 0.0218 yields sin 2 # 12 — 0.341, sin 2 #23 ~ 0.395 or 0.605 which 
are within the experimentally preferred 3a ranges although the solar mixing angle is 
somewhat close to the present 3a upper limit [4]. 


The groups of type C can be generated by three generators a, c, and d, and the generator 
b is dropped. Consequently the diagonalization matrices V) and V v can be A"i and A 3 
shown in Eqs. (3.13, 3.22). As a result, the lepton flavor mixing matrix would be a unit 
matrix for (Vi,V u ) = (A/"i,A/"i), democratic mixing for (V),V/) = (A), A 3 ), (JV3, A/i) or a 
form of Eq. (3.45) for (Vi,V u ) = (A 3 , A 3 ). As discussed above, all three mixing patterns 
are excluded by neutrino oscillation data. Therefore if lepton flavor mixing is completely 
determined by residual symmetries that are subgroups of flavor symmetry, type C groups 
such as the typical T n and A(3n 2 ) series are not suitable as a flavor symmetry. However, in 
the indirect approach [8] where flavor symmetry is completely broken, viable lepton mixing 
may be obtained from a flavor symmetry group of type C for some special vacuum alignments. 


3.3 Predictions of lepton mixing matrix for Majorana neutrino 

If neutrinos are assumed to be Majorana particles, the residual subgroup Q v should be a 
Klein group Z 2 x Z 2 , while Qj can still be any Abelian subgroups of Qjr with order not less 
than 3. For a type D flavor symmetry group, its Z 2 elements are given by 

abc Xl d y \ a 2 bc X2 d y \ bd y \ (3.53) 

where the integers 2 + yi (i = 1 , 2 ) satisfy 

kx\ — ryi = 0 (mod m), {k + l)x 2 — ry 2 = 0 (mod m). (3.54) 

There are another three order 2 elements 

d n/2 , c m/2 d kn/2 , c m/2 d {k+1)n/2 (3.55) 
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Figure 2: Correlations between the mixing angles for the viable mixing pattern in Eq. (3.47), where the 
flavor symmetry is a group of type D. 


for even n. Then m should be even as well since r = m/n = 1 or 3. It is wellknown that a 
Klein group is generated by two order 2 elements. We find that groups of type D have the 
following four Klein subgroups if n is divisible by 2. 

iff = {1 ,d n/2 ,c m/2 d kn/2 ,c m/2 d {k+1)n/2 }, 

Kf } = {1 ,d n/2 ,bd y3 ,bd V3+n/2 }, 

ilf = {1 ,c m/2 d kn/2 ,abc xi d yi ,abc xi+m/2 d yi+kn/2 }, 

Af = {1, c m / 2 d {k+1)n/2 , a 2 bc X2 d y \ a % c ^+m/ 2 d y 2 +{k+i) n / 2 } _ ( 3 . 56 ) 

Notice that k { 1J is a normal subgroup while the latter three Klein subgroups are conjugate 
to each other as 

Af = a 2 K^a , Af = aAf a 2 . (3.57) 

If the full Klein symmetry is preserved by the neutrino mass matrix, the unitary transfor¬ 
mation V v that diagonalizes the neutrino mass matrix can easily be obtained by evaluating 
the eigenvalues and eigenvectors of the representation matrices. The explicit form of V v for 
each residual Klein symmetry is given by 


Af ) : V v 

Af } : K 
Af : K 
Af : K 


1 0 0\ 
010 
0 0 1 / 


1 

71 

1 

71 


72 


/ 0 0 72\ 

rfvz - v ry 3 0 ^ Af 2 {s = 0 ,t = y 3 ), 

\ 1 1 0 / 

/ 0 r]~ Xl — rj~ Xl \ 

I 72 0 0 A/" 4 (s = xi,t = yi), 

V 0 1 1 ) 

/ Q ^-X 2 _ ^-X 2 \ 

I 0 1 1 1 -^ 7V 5 (s = X 2l t = y 2 ), 

\72 0 0 / 


(3.58) 


where the constraints in Eq. (3.54) are taken into account. From the general analysis of 
Sec 3.2, we see that the PMNS would be a unit matrix, a democratic mixing matrix, or 
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a rotation matrix in the (23)-plane if the neutrino residual symmetry Q v = K/\ The 
predictions in agreement with the experimental data can be achieved for Qi = (ac s d*) and 
Q v = K^\ K/\ or K^\ and the PMNS matrix is of the trimaximal form 


tfpMNS — 


/ y/2 COS 9 

I — \/2cos ( 9 — |) 
\ —-\/2 cos ( 9 + |) 


1 — y/2 sin 9 \ 

1 \/2sin (0 - |) I , 
1 \/2sin ($ + §)/ 


with the angle 


9 


(k+Vs-rlt-n) g _ k@) 

s„ = /if 1 
ft, = /if > 


(3.59) 


(3.60) 


where m = rn and n is even. The corresponding predictions for the mixing angles are given 
in Eq. (3.51). Then we proceed to discuss the case that the flavor symmetry is a finite 
SU(3) subgroup of type C. The Klein subgroup would be uniquely K/\ the charged lepton 
diagonalization matrix Vi can be Af\ or A 3 , and therefore the lepton flavor mixing is a trivially 
identity matrix or the democratic mixing pattern, both of which are not phenomenologically 
viable. 


3.4 Numerical results 

As we have shown in Secs. 3.2 and 3.3, among all the possible mixing patterns derived 
from the type D groups as a flavor symmetry, only the trimaximal mixing can be in agreement 
with experimental data. However, no realistic mixing patterns can be obtained from the type 
C groups unless the remnant symmetries Qi and Q v of the lepton mass matrices are partially 
contained in the flavor symmetry group or are purely accidental. The explicit form of the 
PMNS matrix Hpmns is presented in Eq. (3.47), and the predictions for the mixing angles 
are shown in Eq. (3.51). Obviously both Hpmns and mixing angles periodically depend on a 
single real parameter 9 with period 2n. Without loss of generality, the fundamental interval 
of the parameter 9 is chosen to be (— 7 r, 7 r ] in the following. 9 can take a series of discrete 
values, and the expressions of 9 values are given in Eqs. (3.48) and (3.60) for Dirac and 
Majorana neutrinos, respectively. Notice that permutations of rows and columns of the 
PMNS matrix are equivalent to a redefinition of 9 listed in Table 4. 

Depending on the values of the parameters k and r, the groups of type D can be classified 
into two categories: D/% for k — 0, r = 1 and d// n for k = 1, r = 3, as shown in Table 1. 

Since D/ r j n is isomorphic to Z 3 x Dn) n if n cannot be divisible by 3, the representation 
matrices of D./j n can be obtained from those of Dn/ by multiplying 1 , e 27 r V 3 ; Q r e~ 2m / 3 . 
Therefore D/J n and Dn/ give rise to the same sets of PMNS matrices for 3 \ n. It is 
sufficient to focus only on the Dn .11 and D/J, 3n , series. 

a) For groups of Dn/ = A (6 n 2 ) 


We have k — 0 and r = 1 in this case. The parameter 9 is determined by the residual 
symmetries as follows: 


9 

9 


2s-2t-p + 2q 2s-2p + q -2 t + p + q 

7 r, --- 7 T, or --- 7 r, tor Dirac neutrinos, 


2 n 

s~t + y 3 
n 


2 n 


- 7 T, 


Xi 


-t + x 2 


n 


n 


2 n 

n, for Majorana neutrinos , 


(3.61) 
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where s, t,p, q, x\, X 2 , 7/3 — 0,1,..., n — 1. Obviously 9 can take the values 

12 1 

6 (mod 27 t) = — 7 r H- 7 r, — 7r H- 7r,..., 7 r-7r, 7 t, for Dirac neutrinos, 

2 n 2 n 2n 

12 1 

9 (mod 27r) = — 7 r H—7r, — 7 r H—7r,..., 7 r-7r, 7 t, for Majorana neutrinos, (3.62) 

7i n n 

where n has to be an even number in the case of Majorana neutrinos. We see that 9 can 
take 4n and 2 n discrete values for Dirac and Majorana neutrinos, respectively. In other 
words, only half of the 9 values allowed for Dirac neutrinos can be achievable for Majorana 
neutrinos. We perform a scan from n — 2 to n — 100, and the possible predictions for the 
mixing angles sin# 13 , sin 2 9 12 , and sin 2 9 2 3 for each ni°i group are plotted in Fig. 3, where we 
require all three mixing angles to be within the 3a intervals [ ]. To clearly see what we can 
obtain from flavor symmetry groups of relative small order, the subfigures with n ranging 
from 2 to 20 are picked out. Here n = 1 is excluded because Dfl = S3 does not have three 
dimensional irreducible representations. We indicate the results for Majorana neutrinos 
with red points. Note that the same results can also be achieved for Dirac neutrinos; i.e., 
there are always black points that overlap with red points. From Fig. 3, we see that only 
five groups of order less than 2000 can give rise to mixing angles that are in agreement with 
experimental data at the 3a level. Since the five points are located at n = 9,11,14,17,18, the 
corresponding flavor symmetry groups are A(486), A(726), A(1176), A(1734), and A(1944), 
respectively. The value of 9 should be fixed at ±j|, ±§§, ±§f, and ±j|, respectively 
to obtain viable mixing angles. Furthermore, we notice that only the last group A (1944) 
allows neutrinos to be Majorana particles. Although A(1176) has Klein subgroups, it gives 
rise to the lepton mixing angles sin 2 9 13 ~ 0.0148, sin 2 9 \ 2 — 0.338, and sin 2 #23 — 0.414 or 
0.586 for Majorana neutrinos, and the reactor mixing angle $13 is outside the experimentally 
preferred 3a range [ ]. Therefore, neutrinos should be Dirac particles in order to obtain a 
viable lepton mixing pattern from A(1176) flavor symmetry. 

b) For groups of D^, 3n , = (Z 9n > x Z- in ,) x S3 


The parameters are fixed to be k = 1 and r = 3. The analytical expressions for 9 are 
given by 

4s — 6f — 3p + 6 q 2s — 3p + 3q 2s — 6t + 3q 
9 = -——:-7r, - 7 r, or -7r, for Dirac neutrinos, 


9 = 


18tt/ 18n' 18n' 

2s — 3t + 37/3 s — x\ s — 3t + x 2 


9 n' 


-7T 


’ 9n' ’ 


7T, or 


9n' 


n, for Majorana neutrinos, (3.63) 


where s,p , x\, x 2 = 0,1,..., 9n' — 1 and t, q,y 3 = 0,1,... 3 71 ' — 1. Hence 9 can take the 
following discrete values: 

12 1 

9 (mod 27t) = — 7 r H- n, — 7 r H- n ,..., 7 r- n, n, for Dirac neutrinos, 

v ’ 18 n' 18 n' 18 n' 

12 1 

9 (mod 27t) = — 7T + ——-7T, — 7 r + —— 7r,..., 7T — ———7r, 7T, for Majorana neutrinos , (3.64) 

Qt)' Q Q in' 


9 n! 


9 n! 


where n! should be divisible by two as well in the case of Majorana neutrinos. Again only 
half of the 9 values for Dirac neutrinos can survive in the Majorana case. The numerical 
results are plotted in Fig. 4. We see that the mixing angles in accordance with experimental 
data can always be achieved for any value of 71' . Moreover, the first eight D^J, 3n , groups 
with n' — 1, 2,..., 8 lead to the same physical predictions, and accordingly the parameter 9 
is fixed at ±j^. For the first three groups = (Z 9 x Z 3 ) x S 3 , = (Z 18 x Z 6 ) x S 3 and 
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n o ^ 

Figure 3: The viable values of sin 6 * 13 , sin 6 * 12 , and sin" 023 for Dn,n flavor symmetry with respect to the 
group index n, where the predicted mixing angles that fall outside of the 3a limits are not shown. To 
distinguish the predictions for the Dirac and Majorana neutrinos, we present the data in different colors, all 
the red and black points refer to the predictions for Dirac neutrinos while the red points are for Majorana 
neutrinos. The experimentally preferred la and 3cr regions are adapted from Ref. [ ]. 


9 = ( Z27 x Zg) x S3 whose orders are less than 2000, only d [ g 6 = (Z 18 x Z 6 ) x S3 comprises 
Klein subgroups, and neutrinos could be Majorana particles. The same values of the mixing 
angles as those of Refs. [17,38] are obtained. I 11 contrast with A(486) for Qjr = A(6n 2 ), 
0 $ = (Z» x Z 3 ) x S 3 with order 162 is the smallest group that can accommodate the 
experimental data. 

3.5 Predictions of quark mixing 

Analogous to what we have done for the lepton sector, the three generations of left- 
handed quarks are assigned to a faithful triplet of the flavor symmetry group Q^, which is 
broken down to remnant subgroups Qu and Qd in the up quark and down quark sectors, 
respectively. Both Qu and Qo can be any Abelian subgroups that can distinguish among the 
three generations. Therefore the possible forms of the lepton mixing matrix worked out in 
Sec. 3.2 for Dirac neutrinos are also applicable to the quark flavor mixing. It is well known 
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Figure 4: The viable values of sin(?i 3 , sin 2 # 12 , and sin 2 0 23 for Dg^, 3n , — (Zg n i X Z^n') x S 3 flavor symmetry 
with respect to the group index ra, where the predicted mixing angles that fall outside of the 3er limits are 
not shown. To distinguish the predictions for the Dirac and Majorana neutrinos, we present the data in 
different colors; all the red and black points refer to the predictions for Dirac neutrinos while the red points 
are for Majorana neutrinos. The experimentally preferred la and 3cr regions are adapted from Ref. [4]. 


that the CKM matrix is predominantly a unit matrix, and the largest mixing is the Cabibbo 
mixing between the first and the second generations [ ]. As a result, a trivial identity CKM 
matrix is a good first order approximation, and it can always be obtained from any flavor 
symmetry group once the residual symmetries Qjj and Qd coincide. It is remarkable that 
a better approximation can be achieved if we choose ( Vjj,Vd ) = (A/2,A/2), (A/4,A/4), or 
(A/ 5 , A/ 5 ), and the corresponding CKM matrix is determined (up to permutations of rows 
and columns) to be of the form 


VcKM 


( cos 9 — sin 9 0\ 
sin 9 cos 9 0 , 

0 0 l) 


(3.65) 
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where 


r (a+1) ' , ^r M ‘^v , w,a6) 

e= J (N t ,Ni) (3.66) 

[ fcM(£=£l±r(ts) Xi W 5 ,W) 

Here we have assumed that the residual subgroups £q/ and £/£> are generated by a a bc s d l and 
a a bc p d q , respectively, or their direct product extensions. Apparently the Cabibbo mixing 
angle can be correctly produced for an appropriate value of 9. The first two rows and 
columns of Eq. (3.65) can be permutated in this framework, and the resulting CKM matrix 
can be obtained by redefining the parameter 9 in Eq. (3.65) as 9 —> — 6 , 9 —y 7t/2 — 9 or 
9 —» 9 — tt/2. If the flavor symmetry group is of type C, Vjj and Vo can be only A/"i or J\f 3 
given in Eqs. (3.13,3.22), the CKM matrix would be a unit matrix, democratic one, or one 
of the form of Eq. (3.45). Only the unit case is viable as a leading order approximation. 
Notice that the Cabibbo mixing of Eq. (3.65) can also be generated from type B flavor 
symmetry groups [42,43], if the first and the second generations of left-handed quark fields 
are embedded into a doublet. Under this assignment, we have considered all the possible 
finite groups with faithful two-dimensional irreducible representation up to order 1000 in 
Sec. 5, and the phenomenologically viable Cabibbo angles that can be achieved are listed 
in Table 9 below. Notice that the type B groups do not have faithful three-dimensional 
irreducible representations. In the following, we shall focus on the type D groups, and the 
series of 9 values can easily be determined for fixed values of r and k. 

a) For groups of Dn} n — A(6 n 2 ) 


The parameters are fixed to be k = 0 and r = 1 in this case, and then 6 can take the forms 


s -2t-p + 2q 
9 = ---7r, 


2s — t — 2p + q 


-n, or 


s + t — p — q 


where s, t,p,q = 0,1, 


2 n 2 n 2 n 

,7i—l. Therefore, the possible values of 9 are 


-7T . 


12 1 

9 (mod 27 t) = —n H-7r, —7r H- n ,..., 7r-7r, 7r. 

2 7 i 2 n 2 n 


(3.67) 


(3.68) 


We show all possible predictions for the CKM matrix element V us = | sin 9\ for each A(6n 2 ) 
group in Fig. 5. The matrix element V us has been quite precisely measured with \V US \ = 
0.2253±0.0008 [l], such that the currently allowed region appears to be a narrow green line. 
If a A(6n 2 ) flavor symmetry group cannot give rise to lepton mixing angles compatible with 
experimental data, the corresponding predictions for |Ws| are indicated with red points. For 
n — 7 or n — 14, we can predict the Cabibbo mixing angle to fulfill |V^ S | = sin ^ ~ 0.2225, 
which is close to the experimental result. Furthermore, the A(6n 2 ) group with n = 14 can 
generate lepton mixing angles within the experimentally preferred 3a regions while it cannot 
for n — 7. 


b) For groups of D^j, 3 , ^ (Z 9n , x Z M ) x S 3 


We have k — 1 and r = 3 in this case, and the expressions for 9 are 


9 = 


s — 2t — p + 2q 


-7T, 


t — p + q 


- 7 r, or 


t - q 


-7T, 


6 n' ’ 6 n' 3 6 n' 3 ^ 

where s,p = 0,1,..., 6 n' — 1 and t,q = 0,1,..., 3 n! — 1. The series of the 6 values are given 
by 

12 1 

9 (mod 2 tt) = -t r + — tt, -tt + — tt, ..., tt - — tt, t r . (3.70) 

bn' bn' bn' 

The results for |V^ S | = |sind| are displayed in Fig. 6. We see that the D^, 3n , group with 
n' — 7 allows us to predict | V us = sin ~ 0.2225, and meanwhile the observed leptonic 
mixing angles can be accommodated. 
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Figure 5: The possible values of |V US | = |sin0| for D„,h — A(6 n 2 ) flavor symmetry group. The green 
band represents the regions allowed by the experimental data [ 1 ], and the blue band denotes the region 
0.1 < |y us | < 0.3. If the present experimental data of lepton mixing angles cannot be accommodated by a 
A(6 n 2 ) group, the corresponding predictions for \V US \ are plotted with red points. 



Figure 6: The possible values of |V US | = |sin0| for D^, 3n , flavor symmetry group. The green band 
represents the regions allowed by the experimental data [1], and the blue band denotes the region 0.1 < 
\Vus\ < 0.3. 


4 U( 3) subgroup series E(3iV 3 ) and mixing patterns 

The present knowledge on finite subgroups of U( 3) that are not the subgroups of SU(3) 
is partial. Therefore, it is not possible to perform a systematical analysis of mixing patterns 
achievable from all U(3) discrete subgroups as we have done in Sec. 3. In this section, we 
shall focus on the £(37V 3 ) group that is a known series of U( 3) finite subgroups [6]. The 
discrete group £(3 N 3 ) is defined by three Z N symmetries Z N , Z' N , and Z^, which commute 
with each other, and their Z 3 permutations. £(3iV 3 ) is isomorphic to ( Z N x Z' N x Z^) x Z 3 . 
If we denote the generators of Z^,Z' N , and Z^ by a, a', and a", respectively, and the Z 3 
generator by b, then the multiplication rules of the £(3iV 3 ) group are given by 

a N = a' N = a" N = 1, b 3 = 1, 

aa 1 = a'a, aa" = a"a, a'a" = a" a', (4.1) 

b 2 ab = a", b 2 a'b = a, b 2 a"b = a'. 

Consequently all the elements of S(3 N 3 ) can be written as 

g = b k a m a' n a" 1 , (4.2) 

with k — 0,1,2 and m, n, l — 0,..., N— 1. In a specific irreducible faithful three-dimensional 
representation, the generators are presented by the following matrices: 

/ioo\ / 100 \ (e o o\ / 010 \ 

a = I 0 1 0 ) , a' = 0 g 0 , a" = 010 , 6= 001, (4.3) 

\00£/ \0 0 1/ \001/ \ 1 0 0 / 
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where g = e 27r4/W . Next, we proceed to discuss the Abelian subgroups of £(3iV 3 ) as residual 
flavor symmetries and the corresponding diagonalization matrices. 


• Q r = (ba m a' n a" 1 ) 

The representation matrix of the generator is 


/ 0 g n 0 
p 3 {ba m a' n a" 1 ) = 0 0 g m 

\g l 0 0 


with eigenvalues 


^(l+m+n), ^2g\(l+m+n), ^ g\{l+m+n) 


Since the eigenvalues are not degenerate, the corresponding diagonalization matrix is 
uniquely fixed and is given by 

2l+m+n) ^,2 i(—2 l+m+n) ^ g\ (— 2l+m.+n) \ 

g \{-l+2m-n) ug \(-l+2m-n) J2 Q \(-l+2m-n) J . (4.6) 

1 1 1 / 


up to permutations and phases of the column vectors. 
Q r = ( b 2 a m a ,n a" l ) 


/ 0 0 g m 
Ps(b 2 a m a' n a" 1 ) = I g l 0 0 

V 0 g n 0 


The eigenvalues are given by 


^(l+m+n) , u 2g\(l+m+n) ^ n \{l+m+n) 




which are also not degenerate, and the unitary transformation V r is determined to be 


T(— l+2m—n) / .q\{— l+2m—n) ,,2 l+2m—n) 


Af 2 : V r — —p g \{l+m-2n) ^2 g\(l+m-2n) ^ gl(l+m-2n) 

\ 1 1 1 


• g r = ( a m a' n a m ) 


The element a m a' n a" 1 is represented by 


/ g l 0 0 

p 3 (a m a n a" 1 ) = 0 g n 0 

V 0 0 g m 


(4.10) 


It has the eigenvalues 


Q\ Q m ■ 


Obviously the representation matrix is diagonal, and therefore we have 


(4.11) 


/ 1°0 

Ar 3 : V r = 0 1 0 

V o o i 


(4.12) 



for I / n 7 ^ m. The eigenvalues shown in Eq. (4.11) could be partially degenerate such 
that only a mixing vector can be fixed, 


( 1 , 0 , 0 ) T for ?n = n,m^l, 
(0,1, 0) T for 771 = 1,771 ^ 71, 
(0,0, 1) T for 1 = 77,771 ^1. 


(4.13) 


If this is the case, the unitary matrix V r can be fully determined through an extension of the 
residual symmetry Q r = ( a m a' n a" l ) to Q r = ( a m a' n a" 1 ) x Q 2 , where two of three parameters 
m, n, and l are equal, and the generator of Q 2 should be connnutable with a m a ,n a ,n . As a 
consequence, Q 2 can only be a subgroup generated by a p a ,s a' n . If p ^ s ^ t, V r is fixed to be 
A 3 by G- 2 , If any two of p, s, and t are equal, another mixing vector can be determined, and V r 
is a unit matrix as well. Here we assume that mixing vectors associated with ( a m a' n a" l ) and 
G 2 are different; otherwise, the above extension procedure should be continued. However, the 
conclusion of V r being a unit matrix remains. We summarize that the unitary diagonalization 
matrix V r can only be of the form of M\ , 2,3 for E(3A r3 ) flavor symmetry. Assuming neutrinos 
to be Dirac particles, we shall consider all possible combinations of Vi and V v and the 
corresponding predictions for the PMNS matrix in the following. 


• (V), K) = (A/ 3 , A/ 3 ) 

Obviously the PMNS matrix is a trivial identity matrix in this case, 



1 0 0 
0 1 0 
0 0 1 



(4.14) 


which is not viable. 


• {Vi, K) — (Ai, A 3 ), (J\T 3 ,J\fi), (A 2 , A 3 ), (A/ 3 , A 2 ) 

The lepton flavor mixing is found to be the democratic pattern: 



(4.15) 



The PMNS matrix is determined to be of the same form as Eq. (3.45), i.e., 



where 


9 


(l + m — 2n — p + 2 s — t)n 



9 


3N 

(l + m — 2n + p + s — 2t)n 


9 


3N 

(—21 + m + n + p — 2s + t) it 


3N 
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9 


(—21 + m + 7i — p — s + 2t)7T 
3 N 


5 


(m — n — p + s)tt 
N 


for (V h V v ) 


(A/*2, A 2 ) • 


Here we have assumed that the residual symmetries Qi and Q v are generated by the 
elements b k a m a' n a" 1 and b k a p a' s a" t ) respectively, with fc, k' = 1, 2. As shown in Fig. 1, 
the observed values of three leptonic mixing angles cannot be obtained from this mixing 
pattern. Therefore, the £(3iV 3 ) group cannot give rise to a viable lepton flavor mixing 
unless the residual flavor symmetries are partially or completely accidental. If we apply 
the S(3A^ 3 ) flavor symmetry to the quark sector, the CKM matrix would be of the form 
of Eqs. (4.14, 4.15, 4.16) as well, and the Cabibbo angle cannot be generated. 


5 Mixing patterns from a scan of finite discrete groups 

In previous sections, we have performed a comprehensive study of the mixing patterns 
that can be derived from all the SU(3 ) subgroups and the known U( 3) subgroup series 
£(31V 3 ). Although the SU(3) finite subgroups are known very well, the classification of the 
finite subgroups of U( 3) is far from complete. Only a few series such as T,(3N 3 ) are known 
so far. There are still lots of finite groups that are subgroups of U (3) but not subgroups 
of SU( 3). The structures of these groups are not known at present, and therefore the 
possible mixing patterns derived from these group cannot be studied analytically. The task 
of classifying all finite subgroups of U( 3) is a daunting challenge. Instead of doing that, we 
shall perform an exhaustive scan over the discrete groups with order less than 2000 with 
the help of the computer algebra system GAP [25], and the associated flavor mixings for all 
possible residual subgroups of each group are discussed. In the following, we shall focus 
on the groups with faithful three-dimensional irreducible representations to which the three 
left-handed leptons are assigned. 

It is worth mentioning that the groups of orders of 768 and 1536 account for a large 
part of our scan program. There are 1090235 groups of order 768 and 408641062 groups 
of order 1536. The groups of order 768 can be straightforwardly investigated one by one. 
However, it would be rather time consuming to scrutinize through all the groups of or¬ 
der 1536, and the required computational resources are tremendous and beyond our ca¬ 
pability. We notice that all of the groups are regularly arranged in the GAP SmallGroup 
Library [26]. The command SmallGroupsInf ormation( .) can be utilized to get the in¬ 
formation about the arrangement of the groups of a given order. For example, by using 
command SmallGroupsInformation(1536) , we can obtain a list of the group information 
about the groups of order 1536. 

gap>SmallGroupsInformation(1536); 

There are 408641062 groups of order 1536. 

1 - 10494213 are the nilpotent groups. 

10494214 - 408526597 have a normal Sylow 3-subgroup. 

408526598 - 408544625 have a normal Sylow 2-subgroup. 

408544626 - 408641062 have no normal Sylow subgroup. 

This size belongs to layer 8 of the SmallGroups library. 

IdSmallGroup is not available for this size. 

From the above output, we know that the groups of order 1536 are arranged according to the 
properties of their subgroups. The groups of order 768 have the same arrangement pattern, 
which is given by 
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gap>SraallGroupsInformation(768); 

There are 1090235 groups of order 768. 

They are sorted by normal Sylow subgroups. 

1 - 56092 are the nilpotent groups. 

56093 - 1083472 have a normal Sylow 3-subgroup 
with centralizer of index 2. 

1083473 - 1085323 have a normal Sylow 2-subgroup. 

1085324 - 1090235 have no normal Sylow subgroup. 

This size belongs to layer 3 of the SmallGroups library. 

IdSmallGroup is available for this size. 

When we search for the groups of order 768, we find that the groups belonging to the first two 
categories do not admit three-dimensional faithful irreducible representations. We speculate 
that finite nilpotent groups and the groups with normal Sylow 3-subgroup have no three- 
dimensional faithful irreducible representations. This conjecture is checked to hold true for 
many other groups although we are not able to prove it in a strict mathematical way now. 
Under this conjecture, the first 408526597 groups of order 1536 can be skipped, and we only 
need to handle the remaining 114465 groups. We find that only 11 groups of order 1536 
have faithful three-dimensional irreducible representations. A(1536) in the series of A(6n 2 ) 
with n = 16 is known to be a group with faithful three-dimensional representations, it is 
really found in our results, and the group identification number is [1536,408544632], Our 
conjecture is verified again. 

We have studied all the groups with order less than 2000 as flavor symmetry, and the 
groups of order 1536 are treated under our conjecture. Only those groups with faithful 
three-dimensional irreducible representations are considered. Note that the irreducible rep¬ 
resentations given in GAP are not necessarily unitary, we must first convert the nonunitary 
GAP representations into unitary representations. The method of how to construct the cor¬ 
responding unitary representations from a given nonunitary representation is presented in 
the Appendix A. Our method is more efficient than the Gram-Schmidt orthogonalization 
approach in Ref. [ ]. Neutrinos are assumed to be Dirac particles, and then the residual 

subgroups Q v in the neutrino sector and Gi,u,d in the charged fermion sectors can be any 
Abelian subgroups of the postulated flavor symmetry group. We wish to completely fix the 
mixing matrix through breaking of the flavor symmetry to residual symmetries in this work. 
As a consequence, for each group, we find all its Abelian subgroups with order greater than 
two, and the corresponding group structures and generators are extracted. The next step 
is to determine the unitary diagonalization matrices of these subgroups. If a subgroup has 
two or more generators, all the representation matrices of the generators should be diago¬ 
nalized simultaneously. If the eigenvalues of the generators of a subgroup are partially or 
completely degenerate such that the corresponding diagonalization matrix cannot be deter¬ 
mined uniquely (in other words, the residual subgroup cannot distinguish three generations 
of leptons or quarks), such a subgroup would be discarded. Notice that the diagonaliza¬ 
tion matrix is determined up to permutations and phases of their column vectors, since no 
predictions can be made for the fermion masses in the present framework. The first row 
of the diagonalization matrix can always be set to be real and non-negative by redefining 
the overall phases of each column, and the three columns can be arranged such that the 
elements of the first row decrease the column index. A large number of duplicates in the 
unitary diagonalization matrix are removed via this rearrangement procedure. 

Furthermore, we consider all possible residual subgroups for Gi and G v (Gu and Gd hi case 
of quark mixing) for each finite group with three-dimensional irreducible representations, and 
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then the predictions for the PMNS (or CKM) matrix can easily be obtained from the unitary 
diagonalization matrices of the chosen remnant subgroups. Using the freedom to exchange 
the rows and columns for the PMNS matrix, we can order the PMNS matrix in the following 
ways: the (13) element is chosen to be the smallest one, and the (12) element is smaller than 
or equal to the (11) element such that the solar mixing angle is smaller than or equal to 
maximal mixing. Subsequently the flavor mixing angles # 12 , # 13 , #23 and Dirac CP violating 
phase 8cp can be straightforwardly extracted. The exchange of the second and the third 
rows of the PMNS matrix leaves both 9 12 and #13 invariant, while it transforms 6*23 into 
7r/2 — $23 and Sqp to 7r + 8cp■ Notice that all the three mixing angles would keep intact and 
the Dirac CP phase would change from 8cp to — 8cp if we embed the three lepton doublet 
fields into the complex conjugate of the three-dimensional representation. To quantitatively 
measure whether and how well the predicted mixing pattern can accommodate the current 
experimental data, we perform a simple y 2 analysis and the y 2 function is constructed in 
the conventional way, 



and (Jij refer to the best fit values and lcr errors, respectively taken from Ref. [ ]. Note that 
y 2 is very sensitive to the prediction of $ 13 , because the reactor angle $ 13 is quite precisely 
measured. The Dirac phase 8cp is weakly constrained by the present data; therefore, we do 
not include the information of 8cp in the y 2 function in Eq. (5.1). 

The explicit forms of the lepton flavor mixing matrices that can be derived from finite 
groups of order smaller than 2000 have been worked out. Nevertheless, they are too long to 
be included in the present paper. We have made them available at our Website [28], where 
the predictions for mixing parameters, y 2 values, and the associated remnant subgroups are 
presented as well. 


5.1 Results for lepton mixing 

We find that 90 groups with order less than 2000 can generate the mixing angles com¬ 
patible with experimental data up to the 3a level, if neutrinos are Dirac particles. On the 
other hand, only 10 groups can give rise to viable lepton mixing for Majorana neutrinos. To 
our surprise, the PMNS matrix turns out to be of the trimaximal form in all cases, i.e., 


with 


1 / y/2 cos 9 V 1 

Upmns = ~y= ( —V^cos ( 9 U — |) 1 
V3 V-V2 cos ($„ + §) 1 


—\/2sin 9 U \ 
\/2 sin (9 V — |) 
\/2 sin (9 U + f) / 


9 


V 


71 7T 57T 27T 

W’ ± 18’ ± 84’ ± 33‘ 


(5.2) 


(5.3) 


The general expressions of the mixing angles in terms of 9 V are given by Eq. (3.51) with the 
correlations of Eq. (3.52). If we replace 9 V with — 9 U , the solar as well as the reactor mixing 
angles remain invariant, while the atmospheric mixing angle becomes its complement angle 
and the Dirac phase 8cp is shifted by n. The same results are obtained if we exchange the 
second and the third rows of the mixing matrix in Eq. (5.2). For the viable values of 9 V in 
Eq. (5.3), the predictions for the mixing angles and the y 2 value can easily be calculated. 
The numerical results are listed in Table 5. We see that the atmospheric mixing angle $23 is 
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o» 

sin 2 0 i 3 

sin 2 0 i 2 

sin 2 

to 

CO 

y 2 (NO) 

A 2 (IO) 

e u > 0 

e u <0 

e u > 0 

e u < 0 

e u > 0 

e u < 0 

±— 

17 

0.0225 

0.341 

0.393 

0.607 

13.000 

17.456 

33.576 

9.638 

± — 
18 

0.0201 

0.340 

0.399 

0.601 

14.174 

18.807 

34.576 

11.731 

I 07T 

^84 

0.0230 

0.341 

0.392 

0.608 

14.314 

18.724 

34.795 

10.605 

| ‘2tt 

0.0239 

0.341 

0.390 

0.610 

17.520 

21.853 

37.616 

13.076 

±— 

^16 

0.0254 

0.342 

0.387 

0.613 

26.766 

30.947 

45.546 

20.408 

| 57r 
96 

0.0177 

0.339 

0.406 

0.594 

27.034 

31.787 

47.092 

25.495 

I Z7T 

^39 

0.0172 

0.339 

0.407 

0.593 

31.476 

36.245 

51.443 

30.145 

I 07T 

78 

0.0267 

0.342 

0.384 

0.616 

38.449 

42.485 

55.436 

29.808 

± — 
90 

0.0163 

0.339 

0.409 

0.591 

39.606 

44.395 

59.421 

38.608 


Table 5: Mixing angles predicted by finite groups up to order 2000, where we require y 2 is less than 50 
in the case of NO neutrino mass spectrum. For Q v = ±yy, ±y|, ±||, ±|f, all three mixing angles agree 
within 3cr with the experimental data [ ]. For the remaining values Q v = ±fg, ±§f, ±|^, ±||, ±^, both 
solar and atmospheric mixing angles are within the 3a intervals while the reactor mixing angle is slightly 
beyond the allowed range. 


in the first octant and the second octant for positive and negative 9 U , respectively. It is worth 
mentioning that the Dirac CP phase Sep is equal to 0 for 9 U > 0 and 7r for 9 U < 0. Most 
importantly, all mixing patterns with y 2 < 100 lead to a trivial Dirac phase, and nontrivial 
Scp can only arise from patterns that do not accommodate the data well; i.e., the y 2 value 
is larger than 100. 

The finite discrete groups that can give rise to the viable mixing pattern in Eq. (5.2) and 
the corresponding residual subgroups are listed in Table 6 , where the notation [g,j\ denotes 
the 71 h finite group of order g as classified in GAP’s SmallGroups library [26]. The superscript 
* indicates that the corresponding group comprises Klein subgroups such that neutrinos could 
be Majorana particles. The subscripts A and A' imply that the corresponding groups belong 
to the series Dn.l — A (6 n 2 ) and Dg r j, 3n , = (Z 9n / x Z 3n i) x S3, respectively. We see that most 
of these groups have the structures (• • •) x Z 3 x Z 2 or (• • •) x Z 3 x Z 4 , and the orders of the 
groups are divisible by 162 except for the A( 6 n 2 ) groups and the groups that are isomorphic 
to the direct product of A( 6 n 2 ) with cyclic groups. Moreover, we notice that different 
groups may lead to the same PMNS matrix, particularly in the case of 0„ = ±7 t/ 18. For a 
given flavor symmetry group, there are generally many possible combinations of the residual 
symmetries Qi and Q u corresponding to the same Dpmns- Regarding the type D group series 
Dn) n , and D^, 3n , with order less than 2000 , we find the numerical results obtained here with 
GAP are exactly the same as those obtained from analytical study in Sec. 3. This is a very 
convincing check to our results. 

It is remarkable that the smallest groups that can generate the leptonic mixing angles 
within 3cr region of global fits are [162,10], [162,12], and [162,14], and neutrinos are required 
to be Dirac particles. Both [162,10] and [162,12] are U( 3 ) but not 577 ( 3 ) subgroups, and 
[162,14] is exactly the type D group D^l = (Z 9 x Z 3 ) x S3. Although the structure of 
both [162,12] and [162,14] is apparently (Zg x Z 3 ) x S 3 , the definitions of the semidirect 
in the two groups are distinct; therefore, they are not isomorphic. These three groups are 
also found in Ref. [18], where a scan of finite groups up to order 200 has been performed 
under the assumption of Dirac neutrinos. These three groups [162,10], [162,12], and [162,14] 
provide a new possibility for flavor model building. It is interesting to construct a model of 
lepton mixing based on any of these three groups. Such a model is still absent as far as we 
know. It may be better to construct such models in the framework of extra dimension theory 


































rather than supersymmetry as usual, since the tiny Dirac neutrino masses can be naturally 
produced in extra dimension through fermion wave-function localization [! 1, 15]. 

If neutrinos are Majorana particles and the residual symmetry Q v is constrained to be 
a Klein subgroup, then ten groups [648,259], [648,260], [648,266], [1296,1820], [1296,1821], 
[1296,1827], [1944,832], [1944,833], [1944,849], and [1944,2363] can predict acceptable lep¬ 
ton mixing angles, and they lead to the same PMNS matrix that is of the trimaximal form 
with 9 = ±7t/ 18. The groups of order 648 and 1296 are related as follows: 

[1296,1820] 9* Z 2 x [648, 259], [1296,1821] 9* Z 2 x [648, 260], [1296,1827] 9* Z 2 x [648, 266] . (5.4) 

This is exactly the reason why they can lead to the same mixing pattern. Notice that 
only two of the above ten groups are SU{ 3) subgroups, [648, 259] is a group of type D^, 3n , 
with n' = 2, and [1944,849] is a A(6n 2 ) group with n = 18. We see that the results 
of Sec. 3.4 for [648, 259] and [1944, 849] are reproduced by the present numerical analysis. 
The finite discrete groups with order less than 1536 has been considered in Ref. [17] under 
the assumption of Majorana neutrinos, only three groups A(6 x 10 2 ), (Z 18 x Z 6 ) x S 3 , and 
A(6 x 16 2 ) with the identification numbers [600,179], [648,259], and [1536,408544632] are 
found viable, and the reactor mixing angle is predicted to be sin 2 9 13 = 0.0288, 0.0201, 
and 0.0254, respectively. All the mixing patterns claimed in [ ] are confirmed by our 

comprehensive analysis [ ]. However, the values of sin 2 9 13 = 0.0288 and 0.0254 are outside 
the present 3cr interval [ ]. As a result, the groups [600,179] as well as [1536,408544632] are 
omitted in Table 6, and only [648, 259] is left. Comparing our results with those of Ref. [ 7], 
we find five additional groups [648,260], [648,266], [1296,1820], [1296,1821], [1296,1827] 
that can also give the experimentally favored mixing angles, if we concentrate on groups 
with order up to 1536. The reason may be that Ref. [17] only considers the finite subgroups 
of SU(3 ) that cannot be written as direct products with cyclic groups. 

We have also studied all the possible leptonic mixing patterns that can be obtained from 
the “exceptional” finite subgroups £(60) = A 5 , £(60 )xZ 3 , £(168) = PSL{ 2,7), £(168) xZ 3 , 
£(36 x 3), £(72 x 3), £(216 x 3), and £(360 x 3) in our systematic scan. Please see our 
Web page [ ] for details. All the patterns presented in Ref. [19] have been validated, and 

we also obtain other mixing patterns. This is because the requirement that the residual 
flavor symmetries Gi and Q u should generate the whole flavor symmetry group is relaxed in 
our work. We find that only the group £(216 x 3) with identification number [648,532] can 
generate the experimentally acceptable lepton mixing angles. If we take Q v = Z 18 , Gi = Z 3 , 
or Qi = Z 3 x Z 3 , the PMNS matrix would be of the trimaximal form shown in Eq. (3.47) 
with 8 = ± 7 t/ 18 . However, Gi and Q v do not generate the original group £(216 x 3), but only 
a subgroup Dg 1 ] = ( Z 9 x Z 3 ) x S 3 that is denoted by [162,14] in the SmallGroups library. 
The same results are found in Ref. [19]. 

If a certain flavor mixing pattern can be generated from a flavor symmetry group Gt, 
generally the groups containing Gt as a subgroup can also give rise to the same mixing 
pattern. It is interesting to investigate whether the viable mixing pattern predicted by a 
group can also been generated from its subgroups. In other words, we want to discuss whether 
a mixing pattern is intrinsically associated with a group or only one of its subgroups. For 
every group listed in Table 6, we see that there are a large number of possible remnant 
symmetry groups Gu and Gi that lead to the same prediction for 6 U . Each set of {Gu,Gi} 
generates either the parent group or one of its subgroups. In particular, we find 16 groups 
whose residual subgroups cannot give rise to the entire group. These 16 groups and the 
groups generated by remnant symmetries are plotted in Fig. 7. A group and its subgroups 
are linked by arrowed lines, and the arrow points to the parent group. For example, the 
residual subgroups Gu = Z 18 , Gi = Z 3 , or Gi = Z 3 x Z 3 of [486,61], which predict the 
trimaximal lepton mixing with 9 V = ± 7 t/ 18 , only generate the subgroup [162,14] rather 
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than [486, 61] itself. This implies that the group [486, 61] does not give us anything more 
than [162,14] as far as lepton mixing is concerned. As a consequence, [162,14] is preferred 
over [486, 61] as a flavor symmetry group, since they lead to the same lepton mixing pattern, 
and yet the former is less cumbersome and more economical from the view of model building. 
For the remaining 74 groups, the residual symmetries can generate the whole groups as well 
as the subgroups. The correlations between different groups can be shown in a figure similar 
to Fig. 7. The resulting diagram is too complex to be included in the present paper, and it 
is available at our Web site [28]. 

Note that the above predictions for lepton mixing would generally be modified by the 
inclusion of higher dimensional operators. Therefore these mixing patterns that slightly 
deviate from the experimental 3a ranges may still be interesting, since accordance with the 
data could easily be achieved in an explicit model once small subleading corrections are 
taken into account. First, we find that the mixing patterns disfavored at the 3c level with 
X 2 < 50 for NO spectrum are still of the trimaximal form, and the corresponding values of 
9 U are given by 


7r 57T 27r 57T 7r 

A —. A —, A —, A —, A— . 
16 96 39 78 20 


(5.5) 


The predictions for the three lepton mixing angles are displayed in Table 5. The associated 
flavor symmetry groups and the residual symmetries are listed in Table 7. One sees that 
both solar and atmospheric mixing angles are compatible with experimental results while 
the measured value of the reactor mixing angle cannot be accommodated without small 
corrections. 

The second phenomenologically interesting mixing pattern is the transpose of the trimax¬ 
imal mixing matrix. The explicit form of the PMNS matrix is given in Eq. (3.41), where the 
second or the third row of the PMNS matrix is (1,1,1) /\/3. As we have shown in Sec. 3.2, in 
general both # 12 and 9 2 3 have to acquire high order corrections in order to achieve agreement 
with the experimental data in this case. For instance, these three groups [162,10], [162,12], 
and [162,14] can predict the PMNS matrix in Eq. (3.41) with 6 = 7 t/ 9 , which gives rise to 
sin 2 0 i 3 — 0.0201, sin 2 #12 ~ 0.399, sin 2 0 23 ~ 0.340, or sin 2 # 23 ~ 0.660. Obviously, 9 \ 3 agrees 
with the experimental result, and small corrections could bring 612 and 0 23 into the exper¬ 
imentally favored regions. The observed values of d 13 can also be produced by A (6 x ll 2 ), 
A (6 x 14 2 ), and A (6 x 17 2 ) groups, which give 9 = 7tt/66, 37t/ 28, and ll 7 r/ 102 , respectively. 

The third and last interesting mixing pattern takes the following form: 


ffpMNS — 


1 

2 


' A t 9 + 5v^) 

-s [\ (3 - V3) 
W! (3->/3) 


1 

y /2 

1 


0 (9 —5y/3) N 
~ Q(3 + V3) 

0A+7ij J 


(5.6) 


This mixing pattern can be generated from the “exception” finite group £(216 x 3) or the 
group with identification number [648, 531], if the residual symmetries Q v = Z 4 and Gj = Z 18 
are preserved. The mixing angles read 


sin 2 $i 3 ~ 0.0142, sin 2 9 12 — 0.254, sin 2 9 23 = 0.4. 


(5.7) 


If the second and the third rows of the PMNS matrix in Eq. (5.6) are exchanged, the 
atmospheric mixing angle becomes sin 2 # 23 = 0.6. In both cases the atmospheric mixing 
angle is compatible with the experimental data at the 3a level. However, the predicted 
values of 9 \ 2 and 6\ 3 are a bit too small, and this could be reconciled with the experimental 
results in a concrete model after small corrections are included. 
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9 U 

GAP-id 

Group Structure 

Residual Symmetries (Qi Q v ) 

±— 

17 

[1734, 5] a 

7 X Zyj} X Z% X Z2 

(Zq Z34) 


[162,10] 

( Z3 X Z2, X Z2 ) X Z2 X Z 2 

(Zg Z 3 X Zq, Zq) 


[162,12] 

(Zg X Zq ) X ZqX Z 2 

(Zg Z is) 


[162,14] A / 

(Zg x Zq) x Z 3 x Z 2 

(Z 3 , zi 1 Z 18 ) 


[324,13] 

( [Z 3 x Z 3 ) x Z 3 x Z4 x Z 3 

(Zg, Z\ 8 Z 3 X Z\ 2 , Z12) 


[324,15] 

(Zg X Z 3 ) X Z 3 X Z4 

(Zg, Z 18 Z 3 q) 


[324,17] 

(Z 3 X Z2 ) X Z2 X Z 4 X Z2 

(Z 3 x Zq, Z 3 , Z3, Zq Z 3 q) 


[324, 68] 

Z 2 x ((Z 3 x Z 3 x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Z\ 8 Z 2 X Zq, Z 3 X Zq, Zq, Zq) 


[324, 70] 

Z 2 x ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Z\ 8 Z 2 X Z\ 8 , Z 18 ) 


[324, 72] 

Z 2 X ((Zg X Z 3 ) X Z 3 X Z 2 ) 

(Z 3 X Zq, Z 3 , Z 8 , Zq 1 Z 2 X Zis, Z\ 8 ) 


[486, 26] 

(Z 2 j x Z 3 ) x Z 3 x Z 2 

(Z 2 Y Z54) 


[486, 28] 

(Z 2 -j x Z 3 ) x Z 3 x Z 2 

(Z 27 Z 54 ) 


[486, 61] A 

(Zg X Zg) X Z 3 X Z 2 

( Z 3 . Zi I Z 18 ) 


[486,125] 

( Z§ X Z 2 X Z 2 ) X Z 2 X Z 2 

(Z 3 X Zg, Z 3 , Z3, Zg | Z 3 X Zq, Z 3 X Z 18 , 
Zq, z is) 


[648,19] 

( Z2 X Z2) X Z2 X Zg X Z2 

(Zg, Z\ 8 , Z 3 Q Z 3 X Z 2 4, Z 2 4) 


[648, 21] 

(Zg x Z 3 ) x Z 3 x Z 8 

(Zg, Zi 8 , Z 36 Z 72 ) 


[648, 23] 

(Z 3 x Z 3 ) x Z 3 x Z 8 x Z 3 

(Z 3 X Zq, Z 3 X Z\ 2 , Z 3 , Z3, Zq, Z\ 2 Z 72 ) 


[648,128] 

Z4 x ((Z 3 x Z 3 x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Z 18 , Z 3 q Z 2 X Zq, Z 2 X Z 12, Z 3 X Zq, 
Z 3 X Z l2 , Zq X Z12, Zq, Zq, Z\ 2 ) 


[648,134] 

Z4 X ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Z\ 8 , Z 3 q Z 2 x Z\ 8 , Z 2 x Z 3 q, Z\ 8 , Z 3 q) 


[648,140] 

Z4 X ((Zg X Z 3 ) X Z 3 X Z 2 ) 

(z 3 X Zq, Z 3 X Z\ 2 , Z 3 , Zq, Zq, Z\ 2 Z 2 X 
Z\ 8 , Z 2 x Z 3 q, Z\ 8 , Z 3 q) 

± — 
- l 18 

[648,259]* A , 

(Zis x Zq) xZ 3 xZ 2 

(Z 3 , Z 3 \ Z 2 x Zq, Z 2 x Z\ 8 , Z‘2, Z 38 ) 


[648, 260]* 

(z 18 x Zq) xZ 3 xZ 2 

(Zg Z 2 X Zq, Z 2 X Z\ 8 , Z\, Z\ 8 ) 


[648, 266]* 

(Zq x Zq x Z 3 ) x Z 3 x Z 2 

(Zg Z 2 X Zq, Z 3 X Zq, Z\ , Zq, Zq) 


[648, 531] 

Z 3 ■ ((Z 3 x Z 3 ) xQ 8 x Z 3 ) = 
((Z 3 x Z 3 ) x Z 3 x Qs) ■ Z 3 

(Zg Zis) 


[648, 532] 

(Z 3 x Z 3 ) x Z 3 x Q 8 x Z 3 

(Z 3 , zi 1 Z 18 ) 


[648, 533] 

(Z 3 x Z 3 ) xZ 3 xQ 8 xZ 3 

(Zg Z 3 X Zq, Zq) 


[810,21] 

Zq x ((Zq x Z 3 x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Z4Q Z 3 X Zq, Z 3 X Z 30 , Zq, Z 30 ) 


[810,23] 

Zq x ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Z 45 Z 18 , Zgo) 


[810,25] 

Zq x ((Zg X Z 3 ) X Z 3 X Z 2 ) 

(Z 3 X Z 15, Z 3 , Zq, Z 15 Z 18, Zgo) 


[972, 29] 

(Z27 x Z2) x Z2 x Z/± 

(Z 27 , Z\ 54 Z 108 ) 


[972, 31] 

(Z 27 X Z 3 ) X Z 3 X Z4 

(z 27 , 7,54 Z 108 ) 


[972, 64] 

(Zg X Zg ) X Z 3 X Z/± 

(Z 3 X Zq, Z 3 , Zq, Zq Zqq) 


[972,210] 

Z 2 x ((Z 27 x Z 3 ) x Z 3 x Z 2 ) 

(z 27 , 7,54 Z 2 X Z 3 4, Z54) 


[972, 212] 

Z 2 x ((Z 27 x Z 3 ) x Z 3 x Z 2 ) 

(Z 2 Y, Z 3 4 Z 2 X Z 3 4, Z54) 


[972, 245] 

Z 2 X ((Zg X Zg) X Z 3 X Z 2 ) 

(Z 3 X Zq, Zq, Zq, Zq Z 2 X Z 18 , Z 18 ) 


[972, 309] 

(Zg X Z 3 ) X Z 3 X Z4 x Z 3 

(z 3 X Zq, Zq X Zg, Z 3 X Z 38 , Z 3 , Zq, Zq, 
Zg, Z\s Z 3 X Z\ 2 , Z 3 X Zqq, Z\ 2 , Zqq) 


[972, 626] 

Z 2 x ((Zg x Z 3 x Z 3 ) x Z 3 x Z 2 ) 

(Zq X Zq, Zq X Zg, Zq X Z 18 , Zq, Zq, Zq, Zg, 
Zl 8 Z 2 X Zq, Z 2 X Zis, Zq X Zq, Zq X Zqq, 
Zq x Zi 8 , Zq, Zq, Zis) 


[1134,95] 

Z 7 x ((Z 3 x Z 3 x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Zqq Zq X Zq, Zq X Z42, Zq, Z42) 


[1134,97] 

Z 7 x ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Zqq Zis, Z U q) 


[1134,99] 

Zy X ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Zq X Z-21, Zq, Zq, Z 2 i Zis, Z 126) 


Table 6: (continued) 
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GAP-id 

Group Structure 

Residual Symmetries (Qi Q v ) 

1296,35] 

(Z 3 x Z 3 ) x Z 3 x Z 16 x Z 3 

(Zg, z 18 , Z 36 , Z 72 Z 3 x Z 48 , Z 48 ) 

[1296, 37] 

(Zg X Z 3 ) XZ 3 X Z\q 

(Zg, Z\g, Z 3 Q, Z72 Z144) 

[1296,39] 

( Z 3 x Z 3 ) xZ 3 x Z 46 x Z 3 

(Z 3 X Zq, Z 3 X Z\2, Z 3 X Z 2 4, Z 3 , Z3, Zq, 
z 12, Z24 Z144) 

[1296, 274] 

Zg x (( Z 3 x Z 3 x Z 3 ) xZ 3 x Z2) 

(Zg, Z is, Z36, Z72 Z 2 X Zq, Z 2 X Z\ 2 , Z 2 X 
Z24, Z 3 X Zq, Z 3 X Z\2, Z 3 X Z 24 , Zq X Z\ 2 , 
Zq x Z 2 4, Zq, Zg, Z 12 , Z24) 

[1296, 284] 

Z 8 x ((Zg x Z 3 ) xZ 3 x Z 2 ) 

(Zg, Z 18 , Z 3 q, Z 72 Z 2 x Z 18 , Z 2 x Z 3 q, 
Z-2 X Z72, Zis, Z36, Z72) 

[1296, 294] 

Z 3 x ((Zg X Z 3 ) X Z 3 X Z 2 ) 

(Z 3 X Zq, Z 3 X Z12, Z3 X Z24, Z3, Z3, Zq, 
Z 12, Z24 Z 2 x Z is, Z 2 x Z 3 q, Z 2 x Z72, Zis, 
Zqq, Z 72) 

[1296, 688] 

( x Z§ ) x Zq x Z 4 x Zq 

(Z 3 x Zq, Z 3 , Z3, Zq Z 2 x Z 4 , Z 2 x Z 12 , 
Z 2 x Z 3 q, Z 36 ) 

[1296, 689] 

(Zis x Zq) xZ 3 xZ 4 

(Zg, Z 18 Z 2 X Z 4 , Z 2 X Z 12, Z 2 X Z 3 Q, Z 3 q) 

[1296, 699] 

( Zq X Zq ) X Zq X Z 4 X Zq 

(Zg, Z 18 Z 2 x Z 4 , Z 2 x Z 12 , Z 3 x Z 12 , 
Zq X Z 12, Z12) 

[1296,1820]* 

Z 2 X ((Z\g X Zq) X Z 3 X Z 2 ) 

(Z s X Zq, Z 3 , Z 3 , Zq Z 2 X Zq, Z 2 X Z\g, 
Zl x Zq, Zl x Z 18 , z|, Z|, Z 18 ) 

[1296,1821]* 

Z 2 x ((Zis x Zq) x Z 3 x Z 2 ) 

(Zg, Z 18 z 2 x Z 6 , Z 2 x Z 18 , Z| x Z 6 , Z2 x 
Zis, Z\, Z\, Z i 8 ) 

[1296,1827]* 

Z 2 X ((Zq X Zq X Z 3 ) X Z 3 X Z 2 ) 

(Zg, Z 18 Z 2 xZ 6 , Z 2 xZg, Z'^xZq, Z 3 xZq, 
Zl Zl, Zq, Zl) 

[1296, 2893] 

Z -2 x ((Z 3 x Z 3 ) xZ 3 xQ 8 x Z 3 ) 

(Z3 x Zg, Z 3 , Z3, Zg Z2 x Zis, Zis) 

[1296, 2894] 

Z-2 x (Z 3 -((Z 3 x Z 3 ) x Q s x Z 3 ) = 
((Z 3 x Z 3 ) x Z 3 x Q 8 ) ■ Z 3 ) 

(Z9, Zis Z2 x Zis, Zis) 

[1296, 2895] 

Z 2 x ((Z 3 x Z 3 ) x Z 3 x Qs x Z 3 ) 

(Z9, Zis Z2 x Zg, Z 3 x Zg, Zg, Zg) 

[1458,615] 

(Z 8 i x Z 3 ) x Z 3 x Z 2 

(Zgi Z i 62 ) 

[1458,618] 

(Zg 1 x Z 3 ) x Z 3 x Z 2 

(Zgi Zi 62 ) 

[1458, 659] ^ 

(Z 27 X Zg) X Z 3 X Z 2 

(■^ 3 , Z'i Zis, Z54) 

[1458, 663] 

(Z-2Y X Zg) X Z 3 X Z 2 

(Z31 Z 3 Zis, Z 54 ) 

[1458, 666] 

(Z 27 X Zg) X Z 3 X Z 2 

(■^ 3 , -^3 ^ 18 , ZqI) 

[1458,1095] 

(Z 27 x Z 3 x Z 3 ) x Z 3 x Z 2 

(Z3XZ9, Z3XZ27, Z 3 , Z3, Zg, Z27 Z 3 XZq, 
Z 3 x Zis, Z 3 x Z54, Zg, Zis, Z54) 

[1458,1371] 

(Zq X Zq X Zq) X Zq X Z2 

(Z 3 X Zg, Z 3 , Z3, Zg Z3 X Zg, Z3 X Zis, 
Zq, Z is) 

[1620,29] 

Zq X ((Zq X Zq) X Zq X Z 4 X Zq) 

(Zg, Z 18, Z45, Z90 Z 3 X Z12, Z 3 X Zgo, Z12, 
Zqo) 

[1620,31] 

Zq x ((Zg X Z 3 ) X Z 3 X Z A ) 

(Zg, Z 18, Z45, Z90 Z 3 g, Ziso) 

[1620,33] 

Zq X ((Zq X Z%) X Zq X Z 4 X Zq) 

(Z 3 x Zg, Z 3 x Z i5 , Z 3 x Z 30 , Z 3 , Z'I Zq, 
Z15, Z30 Z 36 , Ziso) 

[1620,182] 

Z\q x ((Z 3 x Z 3 x Z 3 ) x Z 3 x Z 2 ) 

(Z9, Zis, Z45, Zgo Z 2 X Zq, Z 2 X Z30, Z3 X 
Zq, Z 3 x Z 30 , Zg x Z 30 , Zg, Zq, Z 30 ) 

[1620,184] 

Z\g x ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Z9, Zis, Z45, Z90 Z 2 X Zis, Z 2 X Zgo, Zi8, 
Zgo) 

[1620,186] 

Z10 x ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Z 3 x Zg, Z 3 x Z i5 , Z 3 x Z 30 , Z 3 , Z'I Zq, 
Z15, Z30 Z2 X Zis, Z 2 X Zgo, Zis, Zgo) 

[1782,21] 

Z\\ x ((Z 3 x Z 3 x Z 3 ) x Z 3 x Z 2 ) 

(Zg, Z99 Z 3 x Zg, Z 3 x Zgg, Zg, Zqq) 
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[1782, 23] 

Z\\ X ((Zg X Z3) X Z 3 x Z 2 ) 

(Zg, Zgg Z is, Z i 98 ) 


[1782, 25] 

Z\\ x ((Zg x Z 3 ) x Z 3 x Z 2 ) 

(Z 3 x Z 33 , Z 3 , Z3, Z 33 Zis, Zigs) 


[1944, 35] 

(Z 2 7 x Z 3 ) xZ 3 xZ 8 

(Z- 27 , Z 3 4 , ZiQ 8 Z 216 ) 


[1944, 37] 

(Z 2 7 x Z 3 ) x Z 3 x Z 8 

(Z 2 7, ^54, ^108 Z 2 iq) 


[1944, 70] 

(Zg X Zg) X Z 3 X Z 8 

(Z 3 X Zq, Z 3 X Z\ 2 , Z 3 , Z3, Zq, Z\ 2 Z 72 ) 


[1944, 362] 

Z4 X ((Z 2 7 x Z 3 ) x Z 3 x Z 2 ) 

(z 27 , Z 54, Z 108 Z 2 X Z 54, Z 2 X Z\q 8 , Z 54, 

^ 108 ) 


[1944, 368] 

Z4 x ((Z 2 7 x Z 3 ) x Z 3 x Z 2 ) 

(z 27 , Z 54, Z\q 8 Z 2 X Z54, Z 2 X Z 10 8, Z 54, 
Zws) 


[1944, 544] 

Z4 X ((Zg X Zg) X Z 3 X Z 2 ) 

(Z 3 X Zq, Z 3 X Zi 2 , Z 3 , Z 3 , Zq, Z \2 Z 2 X 
Z\s, Z 2 x Z 3 q, Z\ 8 , Z 3 Q 


[1944, 707] 

( Zg X Z2, ) X Z2 X Zg X Zg 

(Z 3 x Zq, Z 3 x Zg, Z 3 x Z 12 , Z 3 x Zis, 
Z 3 x Z 36 , Z 3 , Zl, Zq, Zg, Z 12 , Zis, Z 3 Q 

Z 3 x Z 24 , Z 3 x Z 72 , Z 24 , Z 72 ) 


[1944,832]* 

(Z54 x Zq) x Z 3 x Z 2 

(Z 2 7 Z 2 X Zq, Z 2 X Z\ 8 , Z 2 X Z54, Z 2 , Z54) 


[1944,833]* 

(■ Z 54 X Zq) X Z 3 X Z 2 

(Z 27 Z 2 x Zq, Z 2 x Z 18 , Z 2 x Z 54 , Z'2, Z54) 

±— 

- l 18 

[1944,849]; 

(Z\s x Zi 8 ) x Z 3 x Z 2 

(Z 3 , Z 3 Z 2 x Zq, Z 2 x Z\s, Zl Zis) 


[1944,1771] 

Z 4 X ((Zg X Z 3 X Z 3 ) X Z 3 X Z 2 ) 

(Z 3 x Zq, Z 3 x Zg, Z 3 x Zi 2 , Z 3 x Zis, 
Z3 x Z 3 q, Z 3 , Zl Zq, Zg, Zi 2 , Zis, Z 3 Q 
Z 2 x Zq, Z 2 x Zi 2 , Z 2 x Zis, Z 2 x Z 3 q, 
Z 3 x Zq, Z 3 x Zi 2 , Z 3 x Zis, Z 3 x Z 3 q, 
Zq X Zi 2 , Zq X Zis, Zq X Z 3 q, Zq, Z\, Zi 2 , 
Zis, Z 3 q) 


[1944, 2293] 

Zg ■ ((Z 3 x Z 3 ) x Q 8 x Z 3 ) = 
((Z 3 x Z 3 ) X Z 3 X Qs) ■ Zg 

(Z 2 7 Z 54 ) 


[1944, 2294] 

Zg ■ ((Z 3 X Z 3 ) X Qs x Z 3 ) = 
((Z 3 x Z 3 ) x Z 3 x Qs) ■ Zg 

(Z 2 7 Z54) 


[1944, 2363]* 

(Z\s X Zq X Z 3 ) X Z 3 X Z 2 

(Z 3 X Zg, Z 3 , Zl Zg | Z 2 X Zq, Z 2 X Zi 8 , 

Z 3 X Zq, Z 3 X Zis, Zq X Zi 8 , Z\, Zq, Zq, 

Z\s) 


[1944, 3448] 

(Zg X Z 3 ) X Z 3 X Qs x Z 3 

(z 3 X Zg, Z 3 , Z 3 , Zg \ Z 3 X Zq, Z 3 X Zi 8 , 
Zq, Zis) 

1 57 T 
84 

[1176, 243] a 

{Z\^ X Z\/±) X Zg X Z2 

(Z 3 Z4, Z 28 ) 

i 2n 
33 

[726, 5] a 

(Z\\ x Z\\) x Z 3 x Z 2 

(Z 3 Z 22 ) 

[1452,11] 

( Z\\ x Z\\) x Zg x Z/± 

(Z 3 , Zq Z44) 


[1452, 23] 

Z 2 x ((Zn x Zu) x Z 3 x Z 2 ) 

(Z 3 , Zq \ Z 2 x Z 22 , Z 22 ) 


Table 6: Mixing patterns predicted by finite discrete groups up to order 2000, where the mixing angles are 
required to be in the 3cr ranges of the global fit [4]. The superscript * indicates that the corresponding group 
has Klein subgroups and it allows neutrinos to be Majorana particles. The subscripts A and A' denote the 
groups that belong to the Dn} n — A(6 n 2 ) and 3n , = (Zg n > x Z 3n >) x S 3 group series, respectively. The 
notation (Qi \ Q v ) denotes all possible combinations between the subgroups in the Qi list and the subgroups 
in the Q v list. The explicit forms of the residual subgroups can be found at our Web site [28]. 
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e v 

GAP-id 

Group Structure 

Residual Symmetries (Qi Q u ) 


[384, 568] A 

(Zg x Zg) x Zg x Z 2 

(Zg Z 16 ) 


[768,1085335] 

(Zg x Zg ) x Zg x Z4 

(Zg, Z 6 

Z 2 x Zie, Zig) 


[768,1085727] 

Z 2 x (( Zg x Zg) xZgX Z2) 

(Z3, Zg 

Z-2 x Zig, Z\g) 

±— 

- l 16 

[1152,154124] 

Zg x (( Zg x Zg) x Zg x Z 2 ) 

(Z„ zl 

Zig, Z 48 ) 

[1536,408544632]^ 

(Zi6 x Z 3 q) x Zg x Z 2 

(Z3 Z 2 x Z 4 , Z 2 x Zg, Z 2 x Zie, Z 2 , 
Z 4 , Zg, Z 16 ) 


[1536,408544641] 

(Zg x Zg) x Zg x Zg 

(Z3, z 6 , 

Z 12 Z 2 x Zig, Z 4 x Z 16 , Z 16 ) 


[1536,408544715] 

Z 4 x ((Zg x Zg) x Zg x Z 2 ) 

(Z3, Zq, 

Z\ 2 Z 2 X Z\G, Z4 X Z\G, Z\g) 


[1920, 236349] 

Z5 x ((Zg x Zg) x Zg x Z 2 ) 

(Zg, Z\g 

Z\G, Zgo) 

1 57T 
96 

[1536,408544632] A 

(Zie x Z 3 q) x Zg x Z 2 

(Zg Zg 2 ) 

cn|co 

-H 

[1014, 7] A 

(Zig x Z\g) x Zg x Z 2 

(Zg Z 2 q) 

| 57r 
78 

[1014, 7] A 

(Zi 3 x Z 13 ) x Zg x Z 2 

(Zg z 26 ) 


[600,179] A 

(Z 10 x Z 10 ) x Zg x Z 2 

(Zg Z 4 , 

Z 20 ) 

± n 

20 

[1200,682] 

(Z 10 x Z 10 ) xZ 3 xZ 4 

(Zg, Zq 

Z 2 x Z 4 , Z 2 x Z 2 o, Z 4 , Z 2 o) 

[1200,1011] 

Z 2 x ((Z\qX Zio) xZgxZ 2 ) 

(Zg, z 6 

Z 2 x Z 4 , Z 2 x Z 2 o, Z4, Z 2 q) 


[1800, 687] 

ZgX ((Z w xZ 10 ) xZgxZ 2 ) 

(Z 3 , Zl 

Zg, Zg 2 , Z 2 q, Z 6 o) 


Table 7: The mixing patterns with y 2 less than 50 for NO spectrum generated by finite discrete groups 
up to order 2000. The patterns within the 3cr ranges of data [4] are not included since they have been 
collected in Table 6. The superscript * indicates that the corresponding group has Klein subgroups and it 
allows neutrinos to be Majorana particles. The subscripts A and A' denote the groups that belong to the 
Dn^ n = A(6 n 2 ) and D^, 3n , = ( Z 9n i x Z 3n ') x S 3 group series, respectively. The notation (Qi \ Q v ) denotes 
all possible combinations between the subgroups in the Qi list and the subgroups in the Q v list. The explicit 
forms of the residual subgroups can be found on our Web site [28]. 



Figure 7: The groups generated by the residual symmetries Q u and Qi in the lepton sector. If a subgroup 
is generated, it is linked to the parent group by an arrowed line and the arrow points to the parent group. 
The blue boxes denote that all the possible residual subgroups {G„,Gi} do not generate the entire group. 
The yellow boxes imply that the group can be generated by at least one set of the remnant symmetries. 
The similar figure including all the 90 groups in Table 6 looks quite complex, and it is available at our Web 
site [28]. 
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Figure 8: The values of the Cabibbo mixing angle predicted by discrete flavor symmetry groups with order 
less than 2000. The blue dots represent the results shown in Eq. (5.10), and the red star refers to the 
experimental best fit value. 


5.2 Results for quark mixing 

The quark mixing matrix has been rather accurately measured. The fit results for the 
magnitudes of all nine CKM matrix elements are [1] 

/ 0.97427 ± 0.00014 0.22536 ± 0.00061 0.00355 ± 0.00015 \ 

|Vckm| = 0.22522 ±0.00061 0.97343 ± 0.00015 0.0414 ± 0.0012 . (5.8) 

\ 0.008861q!ooo32 0.0405ig;gg^ 0.99914 ± 0.00005 / 

In this subsection, we shall investigate whether the observed quark mixing in Eq. (5.8) can 
be derived from a flavor symmetry group. Similar to the lepton sector, the left-handed quark 
doublets are assigned to an irreducible triplet representation of the group. We have scanned 
all the Abelian subgroups of each finite group with order less than 2000. However, we do not 
find any group that can give rise to the correct order of magnitude of the CKM matrix shown 
in Eq. (5.8). Furthermore, we relax the constraints on each CKM matrix element, and we 
only require that 0.1 < \V US \ < 0.3 and \V u b\ < |VA| < |I4s| should be fulfilled. Surprisingly, 
the admissible CKM matrices turn out to uniquely take the following form: 


cos 6„ 


0 


sin 8 q 0 


Vckm = sin 9 q cos 9 q 


0 


where the values of the rotation angle 9 q are 


9 q e 


nTrnnTrnnnTTnnTrTr 3n n 3n 
30’ 28’ 26’ 24’ 22’ 20’ 18’ 17’ 16’ 15’ 14’ 13’ 12’ 34’ TT’ 32 


(5.9) 


(5.10) 


That is to say, the Cabibbo mixing angle can be naturally produced from a finite group. 
This point has already been seen from the analysis of the 577(3) subgroups in Sec. 3.5. 
The Cabibbo mixing is a good leading order approximation to the measured value of the 
CKM matrix, as the other two smaller quark mixing angles could arise from higher order 
corrections. For comparison with experimental data, the predicted values of the Cabibbo 
mixing angle are plotted in Fig. 8. It is remarkable that the value 9 q = 7t/ 14 is very 
close to the experimental result. The groups that can generate the discrete set of Cabibbo 
angles in Eq. (5.10) and the associated residual subgroups Gu and Go are summarized in 
Table 8. Notice that interchanging the roles of the subgroups Gu and Gd would transform 
the CKM matrix into its Hermitian conjugate, and therefore the same Cabibbo angle would 
be obtained. Generally the residual symmetries Gu and Gd do not give rise to the entire 
flavor symmetry group. It is interesting that some groups can lead to two values of 9 q 
for distinct residual symmetries. For example, the group A(6 x ll 2 ) with identification 
number [726,5] can predict 9 q = 7 t/22 or 7 r/ll. We see that the group of the smallest order 
is [216,95], it leads to 9 q = 7 t/12 , which deviates a bit from the experimental value, and 
this insignificant discrepancy is expected to be resolved by subleading corrections. What is 
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more, the most favored value 9 q = 7t/14 can be generated from the group A(6 x 7 2 ) with the 
group identification [294, 7]. This result is consistent with our analytical study performed in 
Sec. 3.5. A(6 x 7 2 ) provides a new flavor symmetry for constructing models of quark flavor 
mixing. 

Since the third generation quark is much heavier than the first two generations, it is 
assumed that the first and the second generations of the left-handed quark fields transform 
as a doublet while the third generation is a singlet in the well-known U( 2) flavor symmetry 
theory [46]. The doublet plus singlet assignment for the quark sector is also exploited 
by scanning the finite groups with faithful two-dimensional irreducible representation up 
to order 1000. All the possible remnant subgroups Gu and Go are considered for each 
permissible flavor symmetry group. We find that the observed Cabibbo mixing angle can 
easily be accommodated. The phenomenologically viable value of 9 q , which can be achieved 
in this way, are presented in Table 9 where the constraint 0.223 < \V US \ < 0.227 is imposed. 
The complete list of the admissible values of Q q and the flavor symmetry groups as well as 
the residual symmetries are provided at the Web site [28]. It is remarkable that the dihedral 
group D u with GAP id [14,1] can give rise to the interesting Cabibbo angle 9 q = 7t/14. 


o q 

GAP-id 

Residual Symmetries (Gu Gd) or (Go Gu) 

7T 

30 

[1350,46] A 

(Zio, Z 30 Z 10 , Z 30 ) 

7T 

28 

[1176, 243] a 

(Z2 X Z14, Zf, Z 14 Z4, Z 28 ) 

7T 

26 

[1014, 7] a 

(Z26 Z 2 %) 


[864, 701] A 

(Z 2 x Z 4 , Z 2 x Z 6 , Z 2 x Z i2 , Z|, Z 4 , Zyi Z 8 , Z 24 ) 

7r 

24 

[1728, 2847] 

(Z 2 x Z 4 , Z 2 x Z 42 , Z 4 x Z\2i Z 4 , Z 4 , Z\ 2 Z 2 x Z 8 , Z 2 x Z 24 , 
Z'g , Z 24 ) 


[1728,13052] 

(Z 2 x Z 4 , Z 2 x Zg, Z 2 x Z 12 , Z 2 x Z 4 , Z\ x Zg, Z 2 x Z 42 , Z|, Z|, 
z 4 , Z 12 Z 2 x Z 8 , Z 2 x Z 24 , Z 8 , Z 24 ) 


[726, 5] a 

(Z 22 Z 22 ) 

7T 

22 

[1452,11] 

(Z 44 Z 44 ) 


[1452,23] 

(Z 2 x Z 22 , Z 22 | Z 2 x Z 22 , Z 22 ) 


[600,179] A 

(Z 2 x Zio, Z|, Zio Z 4 , Z 2 o) 

7T 

[1200,682] 

(Z 2 x Z 4 , Z 2 x Z 20 , Z 20 Z 2 x Z 4 , Z 2 x Z 20 , Z 4 , Z 20 ) 

20 

[1200,1011] 

(Z 2 x Z 4 , Z 2 x Z 2 o, Z 4 , Z 2 o Z 2 x Zio, ^2 x Zio, Z'i, Zf, Zio) 


[1800, 687] 

(Z 2 x Z 6 , Z 2 x Zio, Z 2 x Z 30 , Zf, Zio, Z 3 o Z 4 , Z i2 , Z 20 , Z 60 ) 


[486,61] A 

(■^18 ■Z’is) 


[972, 64] 

(-^36 Z 3 q) 


[972, 245] 

(Z 2 x Zi 8 , Zi 8 Z 2 x Zi 8 , Zi 8 ) 


[1458, 659] A ' 

(Z\8, Z 54: Zi 8 , Z 54 ) 


[1458, 663] 

('Z’lS) -^54 ^18) ^ 54 ) 

7T 

18 

[1458, 666] 

(Zis, z 54 Zi 8 , Z 54 ) 


[1458,1371] 

(Z 3 X Zg, Z 3 X Zi 8 , Zg, Z i8 Z 3 X Zg, Z 3 X Z i8 , Zg, Zi 8 ) 


[1944, 70] 

(Z 72 Z7 2 ) 


[1944, 544] 

(Z 2 X Zi 8 , Z 2 X Z 3 g, Zi 8 , Z 3 g Z 2 X Zi 8 , Z 2 X Z 3 g, Z 38 , Z 3 g) 


[1944,849]; 

(Z 2 x Zg, Z 2 x Z i8 , Z'2, Zi 8 Z 2 x Zg, Z 2 x Zi 8 , Z 2 , Zi 8 ), (Z 4 , 
■^12, ■Z’se ^4, ^12, ^36) 

7r 

17 

[1734, 5] a 

(Z34 z 34 ) 


[384, 568] A 

(Z 2 x Z 4 , Z 2 x Z 8 , Z|, Z 4 , Z 8 Zig) 

7T 

[768,1085335] 

(Z 2 x Z 4 , Z 2 x Z 8 , Z 4 x Z 8 , Z 4 , Z 4 , Z 8 Z 2 x Zig, Zig) 

16 

[768,1085727] 

(Z 2 x Z 4 , Z 2 x Z 8 , Z| x Z 4 , Z 2 x Z 8 , Z 2 , Z 2 , Z 4 , Z 8 Z 2 x Zig, 
^ie) 


Table 8: (continued) 
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GAP-id 


[1152,154124] 

[1536,408544632] A 
1536,408544641] 
[1536,408544715] 

[1920, 236349] 

[1350,46] A 
[294, 7] A 
[588,16] 

[588, 39] 

[882, 38] 

[1176,57] 

[1176,99] 

[1176, 243] a 
[1470,13] 

[1764,102] 
[1764,171] 
[1014, 7] A 
[216,95] A 
[432, 260] 

[432, 538] 

[648, 259];, 

[648,260]* 

[648,266]* 

[648, 563] 

[864, 701] A 

[864, 703] 

[864, 2348] 

[1080, 298] 

[1296.688] 

[1296.689] 
[1296,699] 

[1296.1820] * 

[1296.1821] * 
[1296,1827]* 


Residual Symmetries (Qu Qd) or (Q D Qu ) 


[7i2 x Z4, Z 2 x Zg, Z2 x Zg, Z2 x Z\ 2, Z2 x Z24, Z 2 , Z4, Z$, Z\ 
■^24 | ZiQ, Z 48 ) 

(Z 2 x Z4, Z 2 x Z s , Z 2 x Z 16 , Zf, Z 4l Z 8 , Z 16 \ Z 2 x Z 4 , Z 2 x Z 
Z2 x Z 16 , Zf, Z4, Zg, Z 16 ), (Z 32 | Z32) 

(Z2 x Zg, Z4 x Z§, Z§, Z| j Z2 x Zi6, Z4 x Zi6, Zie) 

(Z2 x Z4, Z2 x Zf, Z2 x Zg, Zf x Z4, Zf x Zg, Z4 x Zg, Z2 x Z4 x Z 
-Zf, ^2) -^4; ^4) ^8 | ^2 X Z 16 , Z 4 X Z 16 , Z 16 ) 

(Z2 x Z4, Z2 x Z§, Z2 x Z10, Z2 x Z20, Z2 x Z40, Zf, Z4, Z§, Zf 

Z 4 p | Zi6, Zgp) _ 

{Zip, Z 3 q [ Zip, Z 30 ) _ 

(Z14 | Z14) _ 

(Z28 1 Z28) _ 

(Z 2 x Z 14 , Z 14 | Z 2 x Z 14 , Z 14 ) _ 

(Z14, Z42 | Z14, Z42) _ 

(Z56 1 Z 56 ) _ 

(Z2 X Z14, Z2 X Z28, Z14, Z28 I Z2 X Z14, Z2 X Z28, Z14, Z28) 
(Z2 X Z14, Zf, Z14 | Z2 X Z14, Zf, Z 14), (Z4, Z28 1 Z4, Z28) 

(Z14, -^70 | ^14; Z 70 ) _ 

(-^28; -^84 | Z28, Zg 4 ) _ 

(Z 2 x Z 14 , Z 2 x Z 42 , Z 14 , Z 42 | Z 2 x Z 14 , Z 2 x Z 42 , Z 14 , Z 42 ) 

(Z26 1 Z 2 g) _ 

(Z 2 x Zg, Zf | Z4, Z12) _ 

(Z 2 x Z4, Z 2 x Z12 | Z 2 x Z4, Z 2 x Z12, Z4, Z12) 

(Z 2 x Z4, Z 2 x Z 12 , Z 4 , Z 12 | Z 2 x Z 6 , Zf x Z 6 , Zf, Zf) 

(Z 2 x Zg, Z 2 x Zis, Zf, Zis | Z4, Z12, Z 3 g) 

(Z 2 x Z 6 , Z 2 x Z 18 , Zf, Z 18 | Z4, Z 12 , Z 36 ) 

(Z 2 x Zg, Z3 x Zg, Zf, Zg, Z| | Z3 x Z12, Z4, Z12) 

(Z 2 x Zg, Z 2 x Zis, Zf | Z4, Z12, Z 3 g) 

(Z 2 x Z4, Z 2 x Zg, Z 2 x Z42, Zf, Z4, Z42 | Z 2 x Z4, Z 2 x Z 

z 2 x Z 12 , Zf, Z4, Z 12 ), (Zg, Z 2 4 1 Zg, Z24) _ 

(Z 2 x Zg, Z 2 x Z 2 4 | Z 2 x Zg, Z 2 x Z 2 4, Zg, Z24) 

(Z 2 x Z4, Z 2 x Z42 | Z 2 x Z4, Z 2 x Z 6 , Z 2 x Z12, Z| x Z4, Z| x Z 
Z| x Z 12 , Z 4 x Z 12 , Z|, Z|, Z 4 , Z|, Z 12 ), (Z 2 x Z 6 , Z| x Z 
z| x Zg, Z| x Zi 2 , Z2, z| I Z 4 x Zi 2 , Z4, Zf, Zi 2 ) _ 


(^Z 2 X Z 4 , Z 2 X Zi 2 | Z 2 X Z 4 , Z 2 X Zg, Z 2 X Zi 2 , Z 2 X Z 4 , z 2 x z, 

Z\ x Z 12 , Z 4 x Z 12 , Z|, Z|, Z 4 , Z|, Z 12 ), (Z 2 x Z 6 , Z| x Z 

Z 2 2 x Zg, Zf x Zi 2 , Zf, Zf | Z4 x Zi 2 , Z4, Zf, Z 12 ) 

(Z 2 x Zg, Z 2 x Zip, Z 2 x Z30, Zf | Z4, Zj 2 , Z 2 q, Zgp) 

(Z 2 x Z4, Z 2 x Z 12 , Z 2 x Z 36 , Z 36 | Z 2 x Z 4 , Z 2 x Z 12 , Z 2 x Z 3 , 

Z4, Zi 2 , Z3g) 




Z4, Zi 2 , Zgg) _ 

(Z 2 X Z4, Z 2 X Zi 2 , Z 2 X Z3g, Z3g j Z 2 x Z4, Z 2 x Zi 2 , Z 2 X Z3, 

Z4, Zi2, Z 8 q) _ 

(Z 2 x Z4, Z 2 x Z\2) Zg x Zi 2 , Zg x Zi 2 , Zi 2 | Z 2 x Z4, Z 2 x Zi 
Z 3 x Z 12 , Zg x Z 12 , Z4, Z 12 ) 

( rv ^ , rv ry . , ry ry ^ , ry ry ry ry \ ry , , ry ry „ , ry 


Zg x Z 12 , Z 6 x Z 12 , Z4, Z 12 ) _ 

(Z 2 x Z4, Z 2 x Z 12 , Z 2 x Z 36 , Z 4 , Z 12 , Z- 
Zf x Zg, Zf x Z 18 , Zf, Zf, Z 18 ) 

(Z 2 X Z4, Z 2 X Zi 2 , Z 2 X Z 3 g, Z4, Zi 2 , Z; 
Zf x Zg, Zf x Z 18 , Zf, Zf, Z 18 ) 

(Z 2 x Z 4 , Z 2 x Z 12 , Z 3 x Z 12 , Zg x Z 12 , Z 4 , 
zf x Zg, Zg x Zg, Zf, Zf, Zg, Zf) 


"36 I Z2 x Zg, Z 2 x Z 4 
"36 | Z2 x Zg, Z 2 x Zi, 







o q 

GAP-id 

Residual Symmetries (Qu Qd) or (Qd Qu) 

7T 

12 

[1296,2113] 

(Z 2 x Z 4 , Z 2 x Z 12 , Z 2 x Z 36 Z 2 x Z 4 , Z 2 x Z 12 , Z 2 x Z 36 , Z 4 , 

Z\2i Z 3 q) 

[1296, 2993] 

(Z 2 x Z 4 , Z 2 x Z\ 2 , Z 2 x Z 3 6 , Z4, Zi 2 , Z 3 q Z 2 x Zq, Z 2 x Z 48 , 
z'i x Z 6 , Zf x Z 18 , Zl zi) 

[1512,500] 

(Z 2 X Zq, Z 2 X Z 14 , Z-2 X Z 42 , Z'2 Z4, Zi2, Z 2 8 , Z 84 ) 

[1728, 2847] 

(Z 2 x Z 4 , Z 2 x Z\ 2l Z4 x Z 42 , Z 4 , Z 4 , Z\ 2 | Z 2 x Z 4 , Z 2 x Z 42 , 
Z4 x Z 12 , Z 4 , Zf, Z 12 ), (Z 2 x Z 8 , Z 2 x Z 24 , Z 8 , Z 2 4 Z 2 x Z$, 
Z 2 x Z 2 4, Z 8 , Z 24 ) 

[1728, 2855] 

(Z 2 x Ziq, Z 2 x Z 48 Z 2 x Z 16 , Z 2 x Z 48 , Z 16 , Z 48 ) 

[1728,13052] 

(Z 2 X Z4, Z 2 X Zq, Z 2 X Z\2, Z' 2 X Z4, Z'2 X Zq, Z2 X Z\2, Z' 2 , Z'2, 
Z4, Z\ 2 Z 2 X Z4, Z 2 X Zq, Z 2 X Z\2, Z\ X Z4, Z\ X Zq, Z\ X Z\2, 
Z 2l Z'2, Z4, Zi 2 ), (Z 2 x Z 8 , Z 2 x Z 2 4, Zg, Z 2 4 Z 2 x Zg, Z 2 x Z 24 , 

Zg, Z 2 4 ) 

[1728,13083] 

(Z 2 x Z4, Z 2 x Zg, Z 2 x Z\ 2 , Z 2 x Z 2 4 | Z 2 x Z4, Z 2 x Zq, Z 2 x Zg, 
Z 2 x Z\ 2 , Z 2 x Z 2 4, Z\ x Z4, Z\ x Zq, Z\ x Zg, Z\ x Z\ 2 , Z\ x Z 24 , 
Z A xZg, Z A xZ l2 , Z A xZ 2 4, Z2, Z'2, Z A , Z\, Zg, Z\ 2 , Z24), (Z 2 xZq, 
Z\ x Z A , Z\ x Zq, Z\ x Zg, Z\ x Z42, Z\ x Z 2 4, Z\, Z\ \ Z 4 x Zg, 
Z4 X Z\2, Z4 X Z-24, Z4, Z\, Zg, Z 12 , Z 2 4) 

[1944, 832]* 

(Z 2 X Zq, Z 2 X Zig, Z 2 X Z 5 4, Z'2, Z 5 4 Z4, Zi 2 , Z 3 Q, Z w g) 

[1944, 833]* 

(Z 2 X Zq, Z 2 X Zig, Z 2 x Z 5 4, Z‘2, Z 3 4 Z4, Z 32 , Z 3 Q, ZiQg) 

[1944,849]; 

(Z 2 X Zq, Z 2 X Z\g, Z'2, Z\g Z4, Z\2, Z 3 q) 

[1944, 2363]* 

(Z 2 X Zq, Z 2 x Z 18 , Z 3 X Zq, Z 3 X Zig, Zq X Zig, Z'2, Zq, Zg, Zig 
z 3 x Zi 2 , Z 3 x Z 36 , Z 4 , Zi 2 , Z 36 ) 

[1944, 2415] 

(Z 2 X Zq, Z 2 X Zig, Z 2 X Z 3 4, Z'2 Z4, Z12, Z 3 Q, Z 108) 

‘dir 

34 

[1734, 5] a 

(^34 -^34) 

7r 

11 

[726, 5] a 

(Z 22 Z 22 ) 

[1452,11] 

(Z44 Z44) 

[1452,23] 

(Z 2 x Z 22 , Z 22 | Z 2 x Z 22 , Z 22 ) 

37T 

32 

[1536,408544632] A 

(Z 2 x Z4, Z 2 x Zg, Z 2 x Ziq, Z‘ 2 , Z A , Zg, Ziq Z 32 ) 


Table 8: The Cabibbo mixing angle predicted by finite discrete groups up to order 2000, where we require 
0.1 < \V US \ < 0.3 and \V u b\ < \V c b\ < |14s|- The three left-handed quark doublets are assumed to transform as 
a faithful irreducible triplet of the flavor symmetry group. The superscript * indicates that the corresponding 
group has Klein subgroups and it allows neutrinos to be Majorana particles. The subscripts A and A' denote 
the groups that belong to the = A( 6 n 2 ) and 3n , = (Z 9n > x Z 3n i ) x S 3 group series, respectively. 
The notation ( Qu \ Qd) denotes all possible combinations between the subgroups in the Qu list and the 
subgroups in the Qd list. The same prediction for the Cabibbo angle would be obtained if the roles of Qu 
and Qd are interchanged. The explicit forms of the residual subgroups can be found on our Web site [28]. 



GAP-id 

o q 

GAP-id 

327T 

447 

[894,3] 

317T 

433 

[866,1] 

3U7T 

419 

[838,1] 

297T 

405 

[810,3] 

287T 

391 

[782,3] 

27tt 

377 

[754,3] 

267T 

[726,3] 

257T 

349 

[698,1] 

247r 

335 

[670,3] 

237T 

321 

[642,3] 

227T 

307 

[614,1] 

217T 

293 

[586,1] 

207T 

279 

[558,5] 

197T 

265 

[530, 3] 

187T 

251 

[502,1] 

7l7T 

990 

[990,11] 


Table 9: (continued) 


38 








































o q 

GAP-id 

o q 

GAP-id 

3577 

488 

[976,6], [976,7], [976,8] 

6977 

962 

[962, 3] 

177r 
237 

[474,5], [948,5], [948,14] 

677r 
934 

[934,1] 

337T 

4^0 

[920,24], [920,26], [920,28] 

t>57T 

906 

[906, 5] 

lo7T 

223 

[446,1], [892,1], [892,3] 

fco7r 
878 

[878,1] 

3l7T 

432 

[864,4], [864,7], [864,8] 

61.77 

850 

[850, 3] 

157T 

209 

[418,3], [836,3], [836,8] 

597r 
822 

[822, 3] 

297T 

404 

[808,4], [808,6], [808,8] 

57tt 

794 

[794,1] 

147T 

195 

[390,11], [780,11], [780,52] 

5577 

766 

[766,1] 

27tt 

376 

[752,5], [752,6], [752,7] 

537r 
738 

[738, 3] 

1377 

181 

[362,1], [724,1], [724,4] 

5177 

710 

[710,5] 

25tt 

348 

[696,25], [696,27], [696,29] 

4977 

682 

[682, 3] 

127T 

167 

[334,1], [668,1], [668,3] 

47tt 

654 

[654, 5] 

357T 

487 

[974,1] 

237r 
320 

[640,6], [640,7], [640,8 


347T 

473 

[946,3] 

457T 

626 

[626,1] 

117T 

153 

[306,3], [612,3], [612,11], [918,3], 
[918,12] 

437T 

598 

[598, 3] 

327T 

445 

[890,3] 

217T 

292 

[584,4], [584,6], [584,8] 


3l7T 

431 

[862,1] 

4l7T 

570 

[570,11] 

10 tt 

139 

[278,1], [556,1], [556,3], [834,4], 
[834,5] 

39tt 

542 

[542,1] 

297T 

403 

[806,3] 

197T 

264 

[528,65], [528,66], [528,67] 

287T 

389 

[778,1] 

377r 
514 

[514,1] 

9n 

125 

[250,1], [500,1], [500,4], [750,2], 
[750,3], [1000,1], [1000,4], [1000,5], 
[1000,6], [1000,8] 

71t r 
986 

[986, 3] 

357T 

486 

[486,1], [972,1], [972,4] 

267T 

361 

[722,1] 

697T 

958 

[958,1] 

17tt 

236 

[472,3], [472,5], [472,7], [944,4], 
[944,5], [944,6], [944,7], [944,9], 
[944, 30] 

677T 

930 

[930,17] 

257T 

347 

[694,1] 

0077 

458 

[458,1], [916,1], [916,4] 

6577 

902 

[902, 3] 

877 

111 

[222,5], [444,5], [444,17], [666,5], 
[666,15], [888,5], [888,37], [888,38], 
[888,39], [888,41] 

637T 

874 

[874, 3] 

3l7T 

430 

[430,3], [860,3], [860,10] 

237r 
319 

[638, 3] 

6 I 7 r 
846 

[846,3] 

157T 

208 

[416,6], [416,7], [416,8], [832,6], 
[832,7], [832,8], [832,64], [832,68], 
[832,531] 

597r 
818 

[818,1] 

227T 

305 

[610,5] 

297T 

402 

[402,5], [804,5], [804,14] 

367T 

499 

[998,1] 

57tt 

790 

[790,3] 

77r 

97 

[194,1], [388,1], [388,4], [582,4], 

[582.5] , [776,1], [776,4], [776,5], 

[776.6] , [776,8], [970,2], [970,3] 

557T 

762 

[762,5] 

347T 

471 

[942,5] 

27tt 

374 

[374,3], [748,3], [748,10] 

2077 

277 

[554,1] 

5377 

734 

[734,1] 

3377 

457 

[914,1] 


Table 9: (continued) 
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o q 

GAP-id 

o q 

GAP-id 

13tt 

180 

[360,25], [360,27], [360,29], 
[720,65], [720,67], [720,68], 
[720,69], [720,71], [720,178] 

32?r 

443 

[886,1] 

5l7T 

706 

[706,1] 

197T 

263 

[526,1] 

257T 

346 

[346,1], [692,1], [692,4] 

317T 

429 

[858,11] 

49-7T 

678 

[678,3] 

67T 

83 

[166.1] , [332,1], [332,3], [498,2], 
[498,3], [664,1], [664,3], [664,4], 
[664,5], [664,7], [830,2], [830,3], 

[996.2] , [996,3], [996,7], [996,9] 

717r 
982 

[982,1] 

47tt 

650 

[650, 3] 

357T 

484 

[968,3], [968,5], [968,7] 

297T 

401 

[802,1] 

23tt 

318 

[318,3], [636,3], [636,11], [954,3], 
[954,7] 

17?r 

235 

[470,3], [940,3], [940,10] 

457T 

622 

[622,1] 

287r 
387 

[774,5] 

67t r 
926 

[926,1] 

117T 

152 

[304,5], [304,6], [304,7], [608,5], 
[608,6], [608,7], [608,11], [608,25], 
[608,99], [912,60], [912,61], 
[912,62], [912,92], [912,93], [912,94] 

657T 

898 

[898,1] 

277r 
373 

[746,1] 

437T 

594 

[594,3] 

lt>7T 

221 

[442,3], [884,3], [884,14] 

21tt 

290 

[290,3], [580,3], [580,14], [870,3], 
[870,7] 

26?r 
359 

[718,1] 

317r 
428 

[856,3], [856,5], [856,7] 

367T 

497 

[994, 5] 

4l7T 

566 

[566,1] 

6l7T 

842 

[842,1] 

57T 

69 

[138,3], [276,3], [276,9], [414,3], 
[414,7], [552,3], [552,23], [552,24], 
[552,25], [552,27], [690,5], [690,7], 
[828,3], [828,9], [828,13], [828,27], 
[966,9], [966,11] 

597r 
814 

[814,3] 

397T 

538 

[538,1] 

347T 

469 

[938, 3] 

297T 

400 

[800,6], [800,7], [800,8] 

247T 

331 

[662,1] 

19tt 

262 

[262,1], [524,1], [524,3], [786,2], 
[786,3] 

33tt 

455 

[910,7] 

147T 

193 

[386,1], [772,1], [772,4] 

377r 
510 

[510,7] 

237T 

317 

[634,1] 

557r 
758 

[758,1] 

32tt 

441 

[882,5] 

97r 
124 

[248,3], [248,5], [248,7], [496,4], 
[496,5], [496,6], [496,7], [496,9], 
[496,30], [744,24], [744,26], 
[744,28], [744,34], [744,36], 
[744,38], [992,4], [992,5], [992,6], 
[992,7], [992,11], [992,18], [992,25], 
[992,96], [992,99] 

717T 

978 

[978,5] 

317T 

427 

[854, 3] 

537T 

730 

[730,3] 

227T 

303 

[606, 3] 

357T 

482 

[482,1], [964,1], [964,4] 

lo7T 

179 

[358,1], [716,1], [716,3] 

697T 

950 

[950,3] 

307T 

413 

[826, 3] 


Table 9: (continued) 
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o q 

GAP-id 

o q 

GAP-id 

17n 
234 

[234,5], [468,5], [468,17], [702,5], 
[702,36], [936,5], [936,37], [936,38], 
[936,39], [936,41] 

21?r 

289 

[578,1] 

67 7T 

922 

[922,1] 

257T 

344 

[688,5], [688,6], [688,7] 

297T 

399 

[798, 23] 

337T 

454 

[454,1], [908,1], [908,3] 

497r 

674 

[674,1] 

657t 

894 

[894, 3] 

47T 

55 

[110,5], [220,5], [220,14], [330,7], 
[330,11], [440,5], [440,34], [440,35], 
[440,36], [440,38], [550,5], [550,13], 
[660,7], [660,11], [660,35], [660,39], 
[770,9], [770,11], [880,5], [880,90], 
[880,91], [880,92], [880,93], 
[880,94], [880,96], [880,199], 
[990,7], [990,11], [990,25] 

637T 

866 

[866,1] 

477r 

646 

[646,3] 

357T 

481 

[962, 3] 

3l7T 

426 

[426,3], [852,3], [852,9] 

27tt 

371 

[742, 3] 

2o7T 

316 

[632,3], [632,5], [632,7] 

&17T 

838 

[838,1] 

197T 

261 

[522,3] 

347T 

467 

[934,1] 

157T 

206 

[206,1], [412,1], [412,3], [618,4], 
[618,5], [824,1], [824,3], [824,4], 
[824,5], [824,7] 

26?r 

357 

[714,11] 

597T 

810 

[810,3] 

117T 

151 

[302,1], [604,1], [604,3], [906,4], 
[906, 5] 

297T 

398 

[398,1], [796,1], [796,3] 

187T 

247 

[494,3], [988,3], [988,10] 

4o7T 

590 

[590,3] 

2577 

343 

[686,1] 

57tt 

782 

[782,3] 




Table 9: The Cabibbo mixing angle predicted by finite discrete groups up to order 1000, where we require 
0.223 < \V US \ < 0.227. The explicit form of the residual subgroups Qu,d can be found at our Web site [ 28 ]. 
Here the first and the second generations of the quark doublet fields are assigned to a faithful doublet of the 
flavor symmetry group, and the third generation transforms as a singlet. 


5.3 Combined results for lepton and quark mixings 

If flavor symmetry acts on both lepton and quark sectors, from Table 6 and 8, we know 
that some finite groups can generate a lepton mixing pattern in the 3a region together with 
an acceptable Cabibbo angle. These interesting groups and the corresponding predictions 
for 9 V and 9 q are presented in Table 10. The associated residual symmetries Q u ,i,u,d can 
be read from Tables 6 and 8 and our Web page [28]. We see that the smallest group is 
A (6 x 9 2 ), which gives rise to the viable trimaximal lepton mixing with 9 V = ±7t/ 18 and 
the Cabibbo mixing angle 9 q = 7t/ 18. The neutrino should be Dirac particles in this case, 
since A(6 x 9 2 ) does not possess Klein subgroups. In addition, the predicted Cabibbo angle 
is a little smaller than the experimental best fit value. However, higher order corrections are 
typically expected to be of order d 2 ~ 0.05 in concrete models [5,8,9]. Such small corrections 
could pull the Cabibbo angle into the experimentally preferred ranges, and the correct order 
of magnitude of the CKM matrix elements V c b and V ts can be achieved. 

The second-smallest groups are [648,259], [648,260], and [648,266], where the first one 
is exactly the type D group D^, 3n , = (Z 9n i x Z 3n i) x S 3 with n' = 2. All these three groups 
lead to 9 V = ±7t/ 18 and 9 q = 7 t/12 . In contrast with the smallest group A(6 x 9 2 ), neutrinos 
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GAP-id 

9u 

h 

GAP-id 

9u 


[486, 61] A 

±— 

^18 

7T 

18 

[1452,11] 

1 27T 

.33 

7V 7T 

22’ 11 

[648,259]* A , 

± — 
- L 18 

7T 

12 

[1452, 23] 

1 27T 

- 31 33 

7T 7T 

22’ 11 

[648, 260]* 

± — 

18 

7T 

12 

[1458, 659] a / 

±— 

18 

7r 

18 

[648, 266]* 

± — 
- L 18 

7T 

12 

[1458, 663] 

± — 
- L 18 

7r 

18 

[726, 5] A 

1 277 

.33 

7T 7T 

22’ 11 

[1458, 666] 

± — 

18 

7T 

18 

[972, 64] 

±— 

18 

7T 

18 

[1458,1371] 

± — 

18 

7r 

18 

[972,245] 

± — 
^18 

7r 

18 

[1734, 5] a 

± — 

17 

7r 

17’ 34 

[1176, 243] A 

1 57T 

84 

7T 7T 

14’ 28 

[1944, 70] 

± — 

18 

7T 

18 

[1296,688] 

±JL 

- L 18 

7T 

12 

[1944, 544] 

± — 

18 

7r 

18 

[1296, 689] 

± — 
-Ms 

7T 

12 

[1944,832]* 

± — 
- L 18 

7T 

12 

[1296, 699] 

± — 
- L 18 

7T 

12 

[1944,833]* 

± — 
- L 18 

7T 

12 

[1296,1820]* 

± — 
- L 18 

7T 

12 

[1944, 849]; 

± — 
- L 18 

7T 7T 

12’ 18 

[1296,1821]* 

± — 

18 

7T 

12 

[1944, 2363]* 

± — 

18 

7T 

12 

[1296,1827]* 

± — 
^18 

7r 

12 





Table 10: Lepton flavor mixing and Cabibbo mixing angle predicted by finite discrete groups with order 
less than 2000, where both left-handed lepton and quark doublets are assumed to transform as irreducible 
three-dimensional representations of the flavor symmetry group. We require the lepton mixing angles to be 
compatible with experimental data at the 3er level and the CKM entries should fulfill 0.1 < \V US \ < 0.3 and 
V u b < | V c b < 114s |- The superscript * indicates that the corresponding group has Klein subgroups and it 
allows neutrinos to be Majorana particles. The subscripts A and A' denote the groups that belong to the 
D ( n}n = A(6n 2 ) and D^, >3n , “ ( Z 9n f x Z 3n >) x S 3 group series, respectively. 

can be either Dirac or Majorana particles in this case, and the Cabibbo angle is predicted 
to be slightly larger than the measured value. This could easily be reconciled with the 
experimental data in models with small corrections. On the other hand, if one would like to 
obtain a Cabibbo angle 9 q = 7t/ 14, which is already consistent with experimental results at 
leading order, then the group A(6 x 14 2 ) with identification number [1176, 243] is a promising 
flavor symmetry. The PMNS matrix would be of trimaximal form with 0 V = ±57t/ 84, and 
all three lepton mixing angles are in the experimentally preferred range for Dirac neutrinos. 

However, subleading corrections to 0 13 are necessary in order to achieve agreement with the 
data if neutrinos are Majorana particles, as shown in Sec. 3.4. 

From Table 10, we see that the values of the parameters 9 V and 9 q are tightly constrained 
and limited, and they can only take the eight possible values: ( 9 u ,9 q ) =(±7t/18, 7t/18), 
(±7t/18,7t/12), (±2tt/33, tt/22), (±2tt/33, tt/11), (±5tt/84, tt/14), (±5tt/84, tt/28), (±tt/17, tt/ 17) 
and (±7 t/17, 37t/ 34). In particular, a multitude of different groups gives rise to the same 
values 9 V and 9 q . For instance, all the groups [486,61], [972,64], [972,245], [1458,659], 

[1458, 663], [1458, 666], [1458,1371], [1944, 70], [1944, 544], and [1944, 849] predict (0„, 9 q ) =(±vr/18, tt/18). 

For each group, generally many distinct residual symmetries Gu, Gi, Gu, and Gd can lead 

to the same prediction for 9 V and 9 q . The group structures generated by every set of 

{Gu, Gi, Gu, Gd} are studied, we find the original group can always be generated, and certain 

subgroups can also be generated in some cases. All the groups generated by the residual 

subgroups {Gu, Gu Gu, Gd} are plotted in Fig. 9. Notice that no new groups besides those 

listed in Table 10 are found. A group and its subgroups are linked by arrowed lines and 

the arrow points to the parent group. As an example, from this figure, we can read that 

the many residual subgroups of [972,245] for (9 u ,9 q ) =(±7t/18, vt/18) give rise to not only 

the whole group [972,245] but also its subgroup [486,61]. This is exactly the reason why 

[972, 245] and [486, 61] lead to the same predictions for 9 V and 9 q . The same values of 9 U and 

9 q obtained from other different groups can easily be understood in this way. From model 
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GAP-id 

e v 

£3 

o q 

$2 

[486,19] 

±— 

- L 18 

[162,10] 

7T 

18 

[18,1], [54,3] 

[486, 21] 

± — 
- L 18 

[162,12] 

7T 

18 

[18,1], [54,3] 

[486, 23] 

± — 

18 

[162,14] A / 

7T 

18 

[18,1], [54,3] 

[648,126] 

± — 
-As 

[162,10], [324,68] 

7 r 

12 

[24,4], [72,26] 

[648,129] 

± — 
-As 

[162,10], [324,68] 

7T 

12 

[24,6], [72,28] 

[648,132] 

± — 
-As 

[162,12], [324,70] 

7T 

12 

[24,4], [72,26] 

[648,135] 

± — 

18 

[162,12], [324,70] 

7 r 

12 

[24,6], [72,28] 

[648,138] 

± — 
- L 18 

[162, 14] a /, [324,72] 

7T 

12 

[24,4], [72,26] 

[648,141] 

± — 
-As 

[162, 14] a /, [324,72] 

7T 

12 

[24,6], [72,28] 

[648,160] 

± — 
- L 18 

[162,10], [324,68] 

7 r 

12 

[24,8], [72,30] 

[648,164] 

± — 

18 

[162,12], [324,70] 

7T 

12 

[24,8], [72,30] 

[648,168] 

± — 

18 

[162,14] A /, [324,72] 

7 r 

12 

[24,8], [72,30] 

[810, 34] 

± — 
- L 18 

[162,10] 

7 r 7 r 

15 ’ 30 

[30,3], [90,7] 

[810,36] 

± — 
- L 18 

[162,12] 

7T 7T 

15 ’ 30 

[30,3], [90,7] 

[810,38] 

± — 
- L 18 

[162,14] A / 

7T 7T 

15 ’ 30 

[30,3], [90,7] 

[972, 22] 

± — 
- L 18 

[162,10], [324,13] 

7T 

18 

[18,1], [36,1], [54,3], [108,6] 

[972, 24] 

± — 
- L 18 

[162,12], [324,15] 

7T 

18 

[18,1], [36,1], [54,3], [108,6] 

[972, 26] 

± — 

18 

[162,14] A /, [324,17] 

7T 

18 

[18,1], [36,1], [54,3], [108,6] 

[972, 203] 

± — 

18 

[162,10], [324,68] 

7T 

18 

[18,1], [36,4], [54,3], [108,23] 

[972, 205] 

± — 
-As 

[162,12], [324,70] 

7 r 

18 

[18,1], [36,4], [54,3], [108,23] 

[972, 207] 

± — 
-As 

[162,14] A /, [324,72] 

7T 

18 

[18,1], [36,4], [54,3], [108,23] 


Table 11: Lepton flavor mixing and Cabibbo mixing angle predicted by finite discrete groups with order 
less than 1000, where the left-handed lepton and quark fields are assigned to triplet and doublet plus 
singlet of the flavor symmetry group, respectively. Q' 3 and Q' 2 refer to the subgroups generated by the three- 
dimensional and two-dimensional representations, respectively, to which lepton and quark fields are assigned. 
We require the lepton mixing angles should be compatible with experimental data at the 3cr level and the 
CKM matrix elements should fulfill 0.1 < \V US \ < 0.3. The subscripts A' denote the groups that belong to 
D 9n',3n' — x Z 3n >) x S 3 group series. 


building perspective, [486,61] is preferred over [972,245], since [486,61] is smaller in size 
than [972, 245] and consequently it is less cumbersome. However, the larger group [972, 245] 
provides more possible symmetry breaking patterns such that it should be more flexible to 
construct models. 

In view of the above appealing model-independent results, it would be interesting to 
construct a model of quark and lepton mixing based on the flavor symmetry A (6 x 9 2 ), 
[648,259], [648,260], [648,266], A(6 x 14 2 ), or other groups listed in Table 10. For instance, 
guided by the remnant symmetries in Tables 6 and 8, we could introduce auxiliary cyclic 
symmetries, some flavon fields, and driving fields in an appropriate way within the context 
of supersymmetry [5,8,9], such that the postulated flavor symmetry is spontaneously broken 
down to the residual subgroups G v .i in the lepton sector and Gu,d hr the quark sector at leading 
order. Then the desired lepton and quark mixing patterns in Table 10 would be produced 
at leading order, and small corrections could appear if subleading order contributions are 
included. 

On the other hand, if the left-handed lepton and quark fields are assigned to triplet and 
doublet plus singlet representations of the flavor symmetry group, respectively, the exper¬ 
imentally acceptable lepton and quark mixing patterns can be produced as well, as shown 
in Table 11. Here only finite groups with order less than 1000 have been exploited. We 
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[1176,243] 


[1296,688] 


[1296,689] 


[1296,699] 


[1452,11| 


[1734,5] 


[1944,70] 


[1944,832] 


[1944,833] 


[726,5] 


[1452,23] 


Figure 9: The groups generated by the residual symmetries Gi, Gu, Qd} for the discrete flavor symmetry 
groups listed in table 10. If a subgroup is generated, it is linked to the parent group by an arrowed line and 
the arrow points to the parent group. 


see that both the three-dimensional and the two-dimensional representations to which the 
lepton and quark fields are assigned only generate some subgroups of the flavor symme¬ 
try group Qjr since both of them are not faithful representations of Qj. It is remarkable 
that the viable trimaximal lepton mixing with 9 U = ±7t/18 and the Cabibbo mixing angle 
9 q = 7t/18, 7t/12, 7t/15 can be produced. The smallest groups [496,19], [496,21], [496,23] 
as well as the second-smallest groups [648,126], [648,129], [648,132], [648,135], [648,138], 
[648,141], [648,160], [648,164], and [648,168] are interesting flavor symmetry groups that 
predict acceptable lepton and quark flavor mixing at leading order. 


6 Conclusions 

Flavor symmetry and its spontaneous breaking is an interesting approach to explain the 
observed lepton or even quark flavor mixing. It has been shown that the mixing patterns 
that can be derived from the discrete flavor symmetry group are rather restricted for neu¬ 
trino being Majorana particles, only the trimaximal mixing pattern can be in agreement 
with experimental data, and Dirac CP phase is trivial [ ]. In addition, it is found that the 

flavor symmetry group should be of order 648 or larger in order to obtain the experimentally 
preferred mixing angles [ I 7,2 ]. Notice that the Majorana neutrino is only a priori assump¬ 
tion; so far we have no sufficient evidence that neutrinos are Majorana particles. In this 
work, we discussed whether new mixing patterns can be obtained and whether the order of 
the required flavor symmetry groups can be degraded under the assumption that neutrinos 
are Dirac particles 

All the finite subgroups of SU(3) have been identified by three mathematicians: Miller, 
Blichfeldt and Dickson [29]. They can be classified into five categories: type A, type B, type 
C, type D, and type E. Since both type A and type B groups do not have an irreducible 
faithful three-dimensional representation, they are not suitable as flavor symmetry groups. 
We perform a systematical and comprehensive analysis of lepton mixing patterns that can 
be obtained from type C and type D groups. We End that type D groups can give rise to 
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lepton mixing patterns compatible with the experimental data at the 3 a level while type C 
groups cannot. The lepton mixing matrix is still found to be of the trimaximal form, and 
the Dirac phase remains trivial. As a consequence, the resulting PMNS matrix shown in 
Eq. (3.47) depends on a single real parameter 6 , and the solar and atmospheric mixing angles 
are predicted to be sin 2 9 12 — 0.341, sin 2 9 23 — 0.395 or 0.605. These decisive predictions can 
be tested by precisely measuring di 2 and 6*23 at future neutrino oscillation experiments such 
as JUNO [ 7], LBNE [ 8 ], and LBNO [49]. Moreover, forthcoming long baseline neutrino 
experiments LBNE [48], LBNO [49], and Hyper-Karniokande [50] can measure the Dirac CP 
phase with a certain precision. If leptonic CP violation is discovered, the present paradigm 
would be excluded. 

The type D subgroup of SU(3) has two independent series, Dn; n — A( 6 n 2 ) and D,^, 3n , = 
( Z 9n f x Z ' in >) x S 3 . The parameter 9 can only take a series of discrete values depending on the 
values of n and n'. The analytical expressions of 6 are presented for both Dirac and Majorana 
neutrinos. Note that Majorana neutrinos require the flavor symmetry group contains a Klein 
group, such that n and n' should be even. We find that only half of the 6 values allowed for 
Dirac neutrinos can be achievable in the case of Majorana neutrinos. Moreover, we present 
the predictions for the mixing angles as a function of n and n'. It is amazing that the D^, 3n , 
group for any n' allows all three mixing angles to be within the experimental 3 cr ranges, and 
the smallest one is Dg 1 ] = ( Z$ x Z 3 ) x S 3 with order 162. In contrast, for the A( 6 n 2 ) group 
series, A (6 x 9 2 ) with order 486 is the smallest group that can accommodate the experimental 
data. 

We further extend the flavor symmetry paradigm to the quark sector. Similar to lepton 
flavor mixing, the CKM matrix arises from the mismatched remnant symmetries Qu and Qd 
to which a flavor symmetry group is broken in the up type and down type quark sectors. It is 
a pity that no SU(3) subgroups can predict the hierarchical mixing pattern among the three 
quarks although the three generations of the left-handed quark doublets are embedded into 
an irreducible three-dimensional representation of the flavor symmetry group. However, we 
can obtain a good leading order mixing pattern where the Cabibbo mixing angle is generated. 
For instance, both groups A (6 x 7 2 ) and D^, 3n , = (Z 9n i x Z 3 n t) x S 3 with n' = 7 can predict 
the Cabibbo angle to fulfill V us = sin ~ 0.2225, which is quite close to the best fit value. 
In addition, the group A (6 x 14 2 ) can generate both acceptable Cabibbo angle \V US \ = sin 
and lepton mixing angles compatible with the data. 

Usually the flavor symmetry is taken to be a subgroup of U (3) in practice, although 
generally it is not necessary. As has been shown in Ref. [ '], the flavor symmetry Qjr can 
possibly be not an U (3) subgroup. However, none of the three-dimensional representations 
of Qjr would be faithful in this case. The same analysis procedures presented in Secs. 3, 4 and 
5 can be performed, and then essentially the mixing patterns originate from the subgroups of 
Qjr that are generated through the triplets. Notice that the groups generated by the triplet 
representation of Qjr are U( 3) subgroups. It is sufficient to take the flavor symmetry group 
to be finite subgroups of f/(3), and no additional mixing patterns would be obtained even 
if we consider the most general flavor symmetry groups of finite order. From the model 
building perspective, at present we still do not know the advantage of flavor symmetry that 
is not [/(3) subgroups. 

Although the structures of the SU(3) subgroups have been clearly solved in mathematics, 
a complete classification of all U{ 3) subgroups is not available so far. As a consequence, 
it seems unfeasible to analytically discuss all the possible flavor mixing patterns arising 
from the U{3) subgroups, as we have done for the SU(3) subgroups. Given the fact that 
a 7/(3) subgroup series E(3 N 3 ) is known, the mixing patterns that can be derived from 
E(3A 3 ) are studied. However, no mixing pattern in the experimentally preferred 3 cr range is 
found if lepton mixing is completely determined by residual symmetries that are subgroups 
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of £(3fV 3 ). Furthermore, we perform a systematical and exhaustive scan over the finite 
discrete groups up to order 2000 with the help of the computer algebra system GAP. Since 
there are a huge number of groups with order 1536, the groups of order 1536 are treated 
under our conjecture that both finite nilpotent groups and the groups with normal Sylow 
3-subgroup have no faithful three-dimensional irreducible representations. This conjecture 
has been checked to hold true for many finite groups, although we cannot prove it now. 
All possible residual subgroups and the corresponding mixing matrices for each group have 
been worked out. In particular, the mixing patterns that can be obtained from the SU(3) 
subgroups of type E are studied. The abundant results of our analysis are available at the 
Web site [28] . We find that the mixing pattern in accordance with experimental data is still 
of the trimaximal form even if neutrinos are Dirac particles, and all mixing patterns with 
X 2 < 100 predict a trivial Dirac phase. A nontrivial 5cp can only arise from patterns that 
do not accommodate the data well. 

We find that 90 groups with order less than 2000 can generate the mixing angles com¬ 
patible with experimental data at the 3a level, if neutrinos are Dirac particles. Nevertheless, 
only 10 groups can give rise to viable lepton mixing for Majorana neutrinos. The smallest 
groups that can predict the leptonic mixing angles within the 3a region are [162,10], [162,12], 
[162,14] and [648,259], [648,260], [648,266] for Dirac and Majorana neutrinos, respectively, 
where [162,14] and [648,259] are SU( 3) subgroups of type D^J, 3ri , with n' = 1,2, and the 
others are U( 3) but not SU(3) subgroups. Therefore the order of the flavor symmetry group 
decreases 4 times by relaxing the condition of Majorana neutrinos, although no new mix¬ 
ing pattern distinct from trimaximal mixing is obtained. These three groups of order 162 
and the associated symmetry breaking patterns provide new starting points for constructing 
neutrino mass models. 

Then we investigate whether it is also possible to derive quark mixing in the same fashion. 
The left-handed quark fields are assumed to transform as a triplet of Qjr. We impose quite 
loose constraints on the CKM matrix elements 0.1 < [V^ s | < 0.3 and \V u b\ < |VA| < \V US \. 
To our surprise, only the Cabibbo mixing between the first two generations of quarks is 
generated and the other two small mixing angles are predicted to be vanishing. We confirm 
that the group A(6 x 7 2 ) can generate the Cabibbo angle 9 q = vr/14, which is close to the 
best fit value, and the smaller group [216, 95] can give 9 q = ir/12. Moreover, the doublet plus 
singlet assignment for the left-handed quark fields is investigated. We perform a scan of all 
finite groups with faithful two-dimensional irreducible representation up to order 1000. The 
experimentally acceptable Cabibbo mixing angles can easily be generated. For instance, the 
dihedral group D i4 can give rise to a Cabibbo angle 9 q = 7t/ 14, which is quite close to the 
measured value. 

Furthermore, we have considered the scenario in which the flavor symmetry acts on 
both lepton and quark sectors. If both the left-handed lepton and quark fields transform 
as an irreducible three-dimensional representation of the flavor symmetry group, we see 
that the smallest group is A (6 x 9 2 ), which gives rise to the phenomenologically viable 
trimaximal lepton mixing with 9 V = ±7t/ 18 and the Cabibbo mixing angle 9 q = 7t/ 18. The 
second-smallest groups are [648,259], [648,260], and [648,266], which predict 9 V = ±7t/ 18 
and 9 q = 7t/ 12. In both cases, the predicted Cabibbo angles slightly deviate from the 
experimental value. However, this insignificant discrepancy is expected to be eliminated by 
higher order corrections in a concrete model. What is more, a larger group A(6 x 14 2 ) can 
produce both experimentally favored lepton mixing angles and acceptable Cabibbo angle 
9 q = 7t/ 14 at leading order. Therefore A(6 x 9 2 ), [648, 259], [648, 260], [648, 266], A(6 x 14 2 ), 
and other groups in Table 10 provide intriguing opportunity for constructing a complete 
model of quark and lepton, where the observed quark and lepton mixing angles can be 
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naturally explained. On the other hand, if the left-handed lepton and quark fields are 
assigned to triplet and doublet plus singlet, respectively, the trimaximal lepton mixing with 
9 V = ±7t/18 together with the Cabibbo angle 9 q = 7t/18 can also be generated by the flavor 
symmetry groups [496,19], [496,21], and [496,23]. 

Scrutinizing all our results [28] carefully, we find that the flavor mixing patterns arising 
from finite discrete groups seem quite restricted and limited even if neutrinos are Dirac 
particles. It is interesting to perform a complete classification of all the mixing patterns 
achievable in this setup, and such an investigation is left for future studies. 


Acknowledgements 

We are grateful to Peng Long for his participation in the early stage of this work. We 
acknowledge Martin Holthausen for helpful discussions about the computer program GAP 
during FLASY2013 in Niigata, Japan. This work is supported by the National Natural 
Science Foundation of China under Grants No. 11275188 and No. 11179007. 


47 


A Converting a non-unitary representation into a uni¬ 
tary representation 

Since GAP employs only rational and cyclotomic numbers, the irreducible representations 
given in GAP may not be unitary. We have to obtain the corresponding unitary representa¬ 
tions before we proceed to determine the mixing matrix from the postulated residual flavor 
symmetries. We can construct the unitary representations from the nonunitary GAP repre¬ 
sentations in the following way. Let p be a linear representation of a finite group Q. First, 
we define __ 

p = J^pHg)p(g)- (A.i) 

g£G 

Obviously P is a Hcrmitian matrix. According to the rearrangement theorem, we have 

p\h)Pp(h) = pHah)p(gh) = P. (A.2) 

geG 

Following the standard procedures in linear algebra, the matrix P can be diagonalized as 

P = UDU - 1 , (A.3) 

where U is a unitary matrix and it is composed of the normalized eigenvectors of P. D 
is diagonal with entries being the eigenvalues of P. Furthermore, we can show that the 
diagonal elements of D fulfill 

D *k = Y Y (Pu(9)) kj (pu(g)) jk = Y Y\ (pu(9)) jk \ 2 > IM 1 )) J 2 = 1, (A.4) 

ffSG j geG j 

where pu(g) = U~ 1 p(g)U, and therefore pu is a representation of Q equivalent to p. The 
matrix P can be decomposed into the square of a Hermitian matrix S, i.e., 

P = S 2 , S = UD 1 ' 2 U~ 1 , S = S^, (A.5) 

where D 1//2 = diag (p/D n , \[D ...). Then we can straightforwardly check that 

p\h) = Sp(h)S~ l (A.6) 

is exactly the desired unitary representation, because it satisfies 

p'\h)p\h) = S~ l (p\h)SSp(h)) S - 1 = S~ 1 S 2 S~ 1 = 1. (A.7) 
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